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Abstract 

This tutorial presents an introduction to the interaction of light and 
matter on the attosecond timescale. Our aim is to detail the theoretical 
description of ultra-short time-delays, and to relate these to the phase of 
extreme ultraviolet (XUV) light pulses and to the asymptotic phase-shifts 
of photoelectron wave packets. Special emphasis is laid on time-delay ex- 
periments, where attosecond XUV pulses are used to photoionize target 
atoms at well-defined times, followed by a probing process in real time 
by a phase-locked, infrared laser field. In this way, the laser field serves 
as a "clock" to monitor the ionization event, but the observable delays 
do not correspond directly to the delay associated with single-photon ion- 
ization. Instead, a significant part of the observed delay originates from 
a measurement induced process, which obscures the single-photon ion- 
ization dynamics. This artifact is traced back to a phase-shift of the 
above-threshold ionization transition matrix element, which we call the 
continuum-continuum phase. It arises due to the laser-stimulated transi- 
tions between Coulomb continuum states. As we shall show here, these 
measurement-induced effects can be separated from the single-photon ion- 
ization process, using analytical expressions of universal character, so that 
eventually the attosecond time-delays in photoionization can be accessed. 
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1 Introduction 

Great advances in experimental physics are being made using attosecond pulses 
of extreme ultraviolet (XUV) radiation, where it is now possible to initiate, 
control and probe electron dynamics in atoms and molecules in real time. A 
natural and perhaps naive question then arises: Is it possible to use attosecond 
pulses to measure the time it takes for a photoelectron to be ionized? It is 
well-known from fundamental principles in quantum theory that time is not a 
direct observable quantity pQ. Here, our aim is to discuss what temporal as- 
pects of photoionization can be measured in state-of-the-art experiments based 
on attosecond XUV pulses and phase-locked infrared (IR) probe fields. This 
tutorial serves as an introduction to the theoretical description of radiative 
processes taking place in extremely short time intervals, with durations compa- 
rable to the temporal response of outer-shell bound electrons within atoms or 
molecules. An interesting point is that the theoretical background provides a 
common framework for XUV pulses that propagate through dispersive media, 
and for non-relativistic photoelectron wave packets that escape from the atomic 
potential. We wish to stress the similarities and differences between these two 
types of wave packets and to discuss the associated light-matter interactions 
occurring on the attosecond timescale. To this end, we shall derive the way 
in which small perturbations in the medium affect the propagation of the light 
pulses and we will detail the present ways of probing photoionization of neutral 
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species in real time. In doing this, special attention must be given to the case 
of long-range Coulomb potentials, which are important for describing electron 
motion in photoionization. The theoretical methods used here span different 
disciplines ranging from ultra-fast optics and non-linear optics, to scattering 
theory, atomic physics and strong-field physics. 

Outline of Tutorial 

The content of this tutorial is divided into seven sections. In Sec. [T] we introduce 
unfamiliar readers to "attophysics" and we present a short overview of the high- 
order harmonic generation process which is used to generate attosecond XUV 
pulses. In Sec. [2] we review the role of the spectral phase for the propagation 
of XUV wave packets. Here, we use the so-called "group delay" to describe the 
coherent superposition of monochromatic waves in the time-domain. Similarly, 
in Sec. [3] we investigate the intrinsic phase of photoelectrons and the correspond- 
ing time-delay of electron wave packets with special attention to the long-range 
Coulomb potential. In Sec. [4] we perform time-dependent perturbation calcula- 
tions for the creation and propagation of photoelectron wave packets generated 
by an attosecond XUV pulse. In Sec. [5] we review the influence of the laser- 
probe field by performing time-dependent perturbation calculations to second 
order in the interaction with the light fields. In Sec. [6] we discuss the state-of- 
the art experimental efforts to measure attosecond time-delays in laser-assisted 
photoionization. Finally, in Section [7] we present our conclusions and outlook. 

1.1 Overview of attosecond physics 

This section provides an "ultra-short" overview of the historical development 
and of the key concepts in attosecond physics. High-order Harmonic Generation 
(HHG), Attosecond Pulse Trains (APT) and Single Attosecond Pulses (SAP) 
will be introduced. We refer the unfamiliar reader to Ref. [21 13 [4j [5] for more 
comprehensive reviews of attosecond physics. 

1.1.1 High-order harmonic generation 

The work on attosecond light pulses stems from a process discovered in the late 
1980s called HHG [HIT]. It was found that a broad plateau in the XUV range, 
containing a comb of almost equally strong, odd harmonics, could be produced 
by focusing an intense ultra-short IR laser pulse into a target of noble gas. A 
sketch of a typical HHG experiment is shown in Fig.[T](a), where first an intense 
IR field interacts with a target of Ar-gas to generate high-order harmonics, then 
these harmonics propagate to a different chamber where they photoionize Ar- 
atoms so that photoelectrons are emitted. 
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Figure 1: A typical HHG experiment: (a) A fraction of the IR laser field is 
converted to XUV through HHG. The XUV field is then filtered out and used 
to photoionize the detection gas. (b) A representative photoelectron spectra 
(full) and XUV photon spectra (dotted) from an HHG experiment using Ar- 
atoms and Al-filtcr. 

With a Titanium:Sapphire laser system (photon energy of Tvuj = 1.55 eV, cor- 
responding to a near-IR laser wavelength of A = 800 nm) using chirped-pulse 
amplification [8], the typical laser intensity used for HHG is II ~ 10 14 W/cm 2 . 
Harmonic conversion efficiencies as high as 10 -6 and photon energies in the soft 
x-ray range have been obtained [HI EH US H21 1131 HI] • Neither the broad comb of 
high-order harmonics, nor their high conversion efficiency can be explained using 
standard theoretical tools of non-linear optics, where the laser field is treated 
as a perturbation |15j . In the perturbative formalism, an increase in harmonic 
order N is always accompanied by a strong reduction in conversion efficiency as 
is evident from the well-known non- linear conversion formula: Iff cx l£ , where 
II <C 1 is the intensity of the laser in a scaled set of variables. As an example, 
the relevant intensity scale for laser-atom interactions is the so-called atomic 
unit of intensity: Z a t = 3.5 x 10 16 W/cm 2 . As mentioned above, optimal con- 
ditions to observe HHG processes are obtained with Titanium: Sapphire laser 
devices operated at II ps 10 14 W/cm 2 , so that, in scaled units Tj, < 10~ 2 . In di- 
rect contrast with such perturbative behaviour, the harmonic comb from HHG 
may contain many harmonic orders of comparable intensity, Jjy ps In> -C II, as 
seen in Fig. [I] (b). 
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A cut-off law for the maximal photon energy in the HHG plateau was found 
numerically to be |16j 



hQmax ~ Ip ~\~ 

where I p is the ionization potential of the atom, and 
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is the so-called ponderomotive energy of the electron in the laser field [T7] . The 
latter energy is associated with the average oscillatory (quiver) motion of a 
free electron driven by a single-frequency field. Using Eq. (JlJ and ^ in the 
limit Up 3> I p , the cut-off increases quadratically for longer laser wavelengths, 
Ti^max oc A 2 , and linearly with the laser intensity, hQ, max oc 1^. This scaling 
has been used to increase the harmonic photon energies into the keV regime 
using laser fields with longer wavelengths (mid-IR wavelengths) jTSJ HH EH HI] • 
Although HHG spectra from simple systems could be reproduced by nu- 
merical simulations based on the resolution of the Time-Dependent Schrodinger 
Equation (TDSE) |16j . a convenient scenario to explain the mechanism of the 
process is based on the so-called three-step model |22j . The three steps are 
illustrated in Fig. [2j 




Figure 2: The three-step model for HHG: (a) When the electric field of the 
laser is large, the electron can (1) tunnel through the Coulomb barrier into the 
continuum. (6) The electron wave packet will then (2) accelerate on almost 
classical trajectories in the laser field. Finally, a fraction of the electron wave 
packet may be driven back to the atom and (3) recombine, emitting a high-order 
harmonic photon. 

A more detailed quantum mechanical theory for HHG was published by Lewen- 
stein and co-workers in 1994 [53], showing that the three-step model could be 



5 



derived from the Schrodinger equation using the Strong Field Approximation 
(SFA). In this way, a more exact cut-off law was derived analytically, 

nn max = i.3i p + 3.2U P , (3) 

corresponding to electrons with quasi-classical (complex) trajectories [2JJ 125] - 
which are closely-related to Feynman's path integral approach to quantum me- 
chanics |26| . This semi-classical model explained why the harmonics were of 
comparable strength, but it also revealed that several semi-classical trajectories 
contribute to the emission of each harmonic. This finding implied that each 
harmonic was emitted in two bursts per half-cycle of the laser, thus leading to 
a complicated temporal structure of the harmonic emission. It has been shown, 
as reviewed in Ref. [27], that the problem of these multiple emission times can 
be circumvented experimentally thanks to the phase-matching properties of the 
macroscopic medium, leading to a "selection" of the emission from the shortest 
family of electron trajectories. 

Also, it was soon understood that such frequency combs of high-order har- 
monics would correspond to light with pulse durations on the attosecond timescale, 
i.e. shorter than any light pulses ever produced, under some restrictive condi- 
tions regarding their relative phases [23 122 • The generation of such "attosecond 
pulses" is a challenging task especially because the experimentalists have to find 
a way of manipulating the relative phase of the harmonics over the large spectral 
bandwidth of the comb [30 . Another major challenge was to devise a scheme 
to measure the duration of the pulses. This was because existing methods used 
in traditional ultra-fast optics could not be directly applied due to the short 
wavelength and relatively low intensity of the harmonics. 

1.1.2 Attosecond pulses 

In 2001, more than a decade after the first HHG process was observed, the 
first Attosecond Pulse Trains (APT) were characterized using a scheme called 
RABITT (Resolution of Attosecond Beating By Interfering Two-photon tran- 
sitions) [31] [32j [33] . As measured in a set of experiments, the attosecond pulse 
duration was ~ 250 as = 250xl0~ 18 s, corresponding to ~ 1/10 of the laser pe- 
riod, and to ~ 1/100 of the laser pulse duration used for HHG. Two attosecond 
pulses were produced per oscillation of the laser field, resulting in a train of ~ 30 
pulses in the total APT. A typical fraction of an APT is depicted in Fig. [3] 
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Figure 3: Attosecond pulse train: The attosecond pulse (shaded grey) duration 
is a fraction of the fundamental laser period (red dashed curve). The pulse 
separation is a half laser period and the sign of the attosecond electric field 
(heavy black curve) changes from pulse to pulse. 

Making the IR pulse shorter leads naturally to fewer attosecond pulses. It 
was soon demonstrated that a Single Attosecond Pulse (SAP) could be gener- 
ated using a few-cycle IR pulse, ~ 5 fs pH]. The basic properties of the SAP can 
be determined using the attosecond-streak camera [33], while more complete 
characterization requires the use of the sophisticated FROG-CRAB scheme 
(Frequency Resolved Optical Gating - Complete Reconstruction of Attosecond 
Bursts) [55]. For such short IR pulses, the Carrier Envelope Phase (CEP) is an 
important parameter that must be controlled. For instance, the generation of a 
SAP can be changed into the generation of two attosecond pulses by carefully 
tuning the CEP of the laser field [37] [35] . One of the most impressive applica- 
tions of SAP was the "direct measurement of light waves" [39] , where the SAP 
was used to map out the oscillating vector potential of a few-cycle optical laser 
pulse. Another useful way of controlling the HHG process is called "polarization 
gating" [301 SU [13] , which relies on the dependence of laser-polarization for 
the yield of HHG. In simple terms, linear polarization implies that the gate is 
open because attosecond pulses are being produced, while elliptical polarization 
suppresses the HHG and closes the gate. 

From a more fundamental standpoint, attosecond pulses are interesting be- 
cause they are shorter than any pulses created by conventional optical lasers. 
Physically, we understand this because the duration of a light pulse can not 
be shorter than its own period. Optical lasers, corresponding to visible and IR 
wavelengths, have a period that is longer than a femtosecond, leading to the 
so-called femtosecond barrier for ultra-short pulses. The power needed to drive 
a laser scales strongly with the photon frequency, P oc w 5 , which effectively pre- 
vents the conventional laser scheme from going beyond optical wavelengths [TH] . 
More formally, a pulse intensity envelope has a minimal duration proportional 
to the inverse of the supporting coherent bandwidth, 
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where Aui (and AA -C A) is the intensity bandwidth of the light. In the case 
of a Gaussian pulse, C = 4 In 2 for the intensity envelope. The minimal pulse 
duration, called the Fourier-limited pulse duration, occurs when all the spectral 
components are compressed, i.e. when the spectral phase is perfectly linear. 
The attosecond pulses produced through HHG can overcome the femtosecond 
barrier not only because they have a much shorter wavelength than optical 
light, but also because they have a broad coherent bandwidth. It is interesting 
to compare the bandwidth of the laser used for HHG with the bandwidth of 
the attosecond pulses being generated. An IR laser pulse used for generation of 
an APT has a typical duration of 30 fs and a bandwidth of 60 meV; while the 
corresponding attosecond pulses are 100 as with a bandwidth of 15 eV. During 
the HHG process the coherent bandwidth is increased by a factor of 250 from 
laser to high-order harmonics. 

Attosecond pulses can be used to ionize atoms and molecules at extremely 
well-defined times, but how do the resulting photoelectron wave packets and 
holes behave after the ionization event? How can the attosecond pulses be 
used to obtain temporal information about physical processes occurring on the 
atomic timescale? Using attosecond pulses to ionize target atoms or molecules 
does not by itself provide any high-resolution temporal information! When the 
photoelectrons are collected in an experiment, we can only determine the value 
of the momentum (or energy) they have gained in the course of the process, 
but not how or when it was acquired. This implies that the temporal ionization 
events must be synchronized and probed with the help of an external "clock" 
that "ticks" itself on the attosecond timescale. The typical period of light in 
the optical regime is on the femtosecond timescale, but one can argue that its 
variation occurs in the sub-femtosecond regime, which then justifies its common 
usage as a probe in attosecond photoionization experiments. In this way, it is 
possible to use a fairly long laser pulse, say tens of femtoseconds, to monitor 
attosecond photoionization events. Clearly, it is important that the probe-field is 
phase-locked with the attosecond pulses, a requirement which is easily satisfied in 
HHG experiments by using fractions of the same laser pulse for both creating the 
attosecond pulses and for probing the photoionization process. Using attosecond 
pulses with phase-locked laser probes, has led to a range of experiments where 
time-delays between photoelectrons from different configurations and systems 
have been measured [HI H3 [46] . Then, the natural question that emerges is 
whether or not such attosecond experiments yield direct access to the delay in 
photoemission, i.e. to the time it takes for an electron to photoionize, or if 
the observed delays should be interpreted in a different way? Before answering 
these questions, we briefly review different pump-probe schemes that have been 
used to experimentally study light -matter interactions in real time. 

1.2 Pump— probe schemes for ultra-fast measurements 

In this tutorial, we provide a theoretical background for the temporal aspects 
of photoionization. The attosecond pulse structure is shorter than the response 
time of any detector or electronic device, but it is possible to obtain temporal 
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information about attosecond photoionization indirectly by investigating coher- 
ent cross-correlation photoelectron spectrograms between the attosecond pulse 
and a weak IR laser probe. In Fig. [4] we present three different kinds of cross- 
correlation techniques. 

(a) 



(b) 




r 



Figure 4: Pump-probe schemes: (a) Traditional pump-probe experiment with 
two pulses separated in time by r. (b) Simultaneous pump-probe experiment 
between a SAP and a few-cycle IR field, (c) Simultaneous pump-probe experi- 
ment between an APT and a monochromatic IR field. The narrow purple area 
represents the attosecond XUV pulse envelope and the broader red area repre- 
sents the one of the probing laser pulse, while the dotted red lines indicate the 
corresponding E-ficld. 

1.2.1 Traditional pump probe experiment 

A typical pump-probe scheme is illustrated in Fig. 0(a). First, a pump pulse 
is used to excite the system at time t = 0. A part of the quantum mechanical 
wave function is then pumped from the initial state to an excited wave packet: 
*W(0) ^ $(i)(o). The syste m can then evolve in a field free environment until 
time t = r, when the probe pulse interacts with the excited wave packet, thus 
changing its state: ^^(r) — > 3K 2 '(t). In experiments, the modification of the 
wave packet leads to the change of the observable quantities as a function of 
the delay between the pump and the probe pulses. Temporal information about 
the field-free propagation of the system, e.g. a molecular vibration, can then be 
extracted in real time by repeating the experiment systematically for different 
delays of the probe pulse. In this way, the intermediate steps in chemical reac- 
tions, the so-called "transition states" , have been investigated in the framework 
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of Zewail's femtochemistry [571 SHI ESI ISO] ■ Clearly, the intuitive interpretation 
of these experiments is limited by the respective pulse durations. Using few-cycle 
optical-laser fields to pump and probe the system is adequate to study nuclear 
motion in molecular systems, typically occurring on the femtosecond timescale. 
Coherent, attosecond XUV pump-probe experiments, where attosecond pulses 
are used for both pumping and probing, hold promise of observing the elec- 
tron dynamics on the attosecond timescale, but they are difficult to implement 
experimentally due to the low probability for absorbing two photons in the 
XUV range, where we assume that at least one photon must be absorbed from 
both the pump and the probe for a meaningful signal. We refer the reader to 
Ref. [5T] (and the references therein) for state-of-the-art experimental efforts on 
XUV-pump and XUV-probe experiments bordering the attosecond timescale. A 
different kind of pump-probe experiments has been carried out using a SAP as 
pump and a few-cycle laser pulse as probe, see for instance Ref. [52| , where the 
Fourier transform of the delay-dependent spectrogram yields quantum beats, 
but also interference between direct and indirect pathways in the ionization. 

1.2.2 Simultaneous pump probe experiment using SAP: "Streak- 
ing" 

A more commonly used attosecond pump-probe configuration is illustrated in 
Fig. [4] (b) , where the pump and probe pulses overlap in time. In this situation it 
is not primarily the field- free system that is of interest, but rather the temporal 
characterization of the pulses or the evolution of the system in the presence of the 
two fields. The pump probe is an attosecond XUV pulse, while the probe pulse 
is a longer (few femtosecond) IR-laser pulse. Strictly speaking, this is a laser- 
assisted photoionization process, where the system is simultaneously pumped 
and probed. Sub-femtosecond temporal information can be gained by repeating 
the experiment at different subcycle delays between the attosecond pump pulse 
and the laser field oscillation of the probe pulse. Clearly, it is essential that 
the two pulses are phase-locked and that the delay can be controlled with sub- 
femtosecond precision. 

Under these conditions, the system can absorb energy simultaneously from 
the two fields i.e. while it absorbs one XUV photon, it is "dressed" by the 
relatively intense IR probe field. Accordingly, a non-perturbative approach is 
needed to account for the effect of the probe. A simplified theory for the in- 
fluence of the probe pulse on the photoelectrons is given by the streak camera 
formalism [3 5) . This simplified interpretation of the photoelectron distribution 
relies on the use of the SFA, which is equivalent to assuming that the photo- 
electron feels the instantaneous laser field while the Coulomb potential of the 
ionic core is neglected. According to the SFA, the streaked electron momentum 
distribution is then shifted as 

Pf( T ) =Po - eA(r), (5) 

where p*o is the unshiftcd, probe-free momentum; e > is the elementary charge 
and A(t) is the probe-field vector potential at the time of ionization. This con- 
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cept of "instantaneous streaking" provides a simple map from time to momen- 
tum of the streaked electrons, which is valid assuming that the electron is free 
shortly after the ionization event, i.e. that it cither escapes from a short-range 
binding potential or is ejected with a high velocity. In practice, the electron 
is not completely free after the ionization event due to the remaining Coulomb 
potential from the ion. Smirnova and co-workers found that the simultaneous 
action of the long-range Coulomb potential and the probing laser field can result 
in a small shift of the streaking process [53], i.e. an uncertainty on the absolute 
delay in the experiment. Until recently [45], these delays were considered too 
small to be accessible in streaking experiments. 

1.2.3 Instantaneous pump probe experiment using APT: "RABITT" 

In Fig. [4] (c), a train of "identical" attosecond pump pulses are probed with 
a monochromatic IR-field. This scheme is often referred to as the RABITT 
method [3"Tl 1521 155] . Provided that the probe pulses are weak and that the fields 
are repeated periodically for many cycles, this scheme can provide equivalent 
information as the streak-camera method described above [Fig. [4] (b)]. We will 
discuss this equivalence in more detail in Sec. [5] in terms of photoelectron wave 
packets, and we suggest that the interested reader should consult Ref. [54] for 
a complete derivation in lowest-order perturbation theory. Higher-order effects 
in the probe field have been considered using the "soft-photon approximation" 
|55j . but the high-level of accuracy that is required for attophysics remains a 
challenge to theory, especially so in the non-perturbative regime. 

The shared periodicity of the APT and the probe field implies that the pho- 
toelectrons will appear on distinct energies corresponding the discrete numbers 
of absorbed photons. In Fig. [5] we display a typical experimental photoelectron 
spectrogram produced by laser-assisted APT ionizing Ar gas. 
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Figure 5: RABITT method: (a) Photoelectron spectrogram over photon energy 
and delay between the APT and the IR field. The offset in the modulation of the 
sidebands contains information about the attosecond pulses and the ionization 
process, (b) Schematic energy diagram over the quantum paths leading to the 
same final energy in sideband (SB) 2q. The experimental data was gathered for 
Ref. [SS]. 

The advantages of the periodic time-structure of APT are rather practical: First, 
it is less demanding experimentally to produce APT and multi-cycle probe- 
fields, than SAP and few-cycle probe fields. Second, the probe field can be 
weaker than for streaking, thereby, causing less side-effects on the system, e.g. 
induced polarization of the core or ionization. Third, the signal is read out 
on zero-background, i.e. in energetic regions where no initial photoelectrons 
appear. Finally, a monochromatic probe field simplifies the analysis because it 
leads to less spectral convolutions in the experimental signal, as we shall discuss 
further in Sec. [U 

In this setup, there is a phase-difference associated with the instantaneous 
probing by the laser field, which can be traced back to the phase-shift of the 
relevant two-photon matrix elements. The probability of the sideband peak 
modulates as 

S 2q = a + (3 cos[2wt - A<j) 2q - A6 2q ] , (6) 

where A(f>2 q = {4>2q+i — 4>2q-i) is the phase difference between the consecutive 
harmonics (2q + l)w and (2q — l)w; and A# 2 q is an intrinsic atomic quantity, 
the so-called "atomic phase" , associated to the difference of the phases of the 
transition amplitudes associated to the distinct quantum paths leading to the 
sideband [32 [33] . The former phase is related to the arrival time of the attosec- 
ond pulses: T0 = A(^2q/(2w), for frequencies f2 ~ 2qco. The atomic phases were 
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studied in detail by Toma and Muller already in 2002, using lowest-order per- 
turbation theory, for photoionization from the 3p state in Argon |57j . Further 
numerical work was performed by Mauritsson, Gaarde and Schafer using single- 
active electron (SAE) effective potentials for Helium, Neon and Argon [35]. The 
two latter numerical calculations were made with the intention of making cali- 
bration curves for attosecond characterization tools such as RABITT. Our aim 
in this tutorial is to discuss the origin of this atomic phase in detail. We will 
identify its fundamental physical components in terms of time-delays and we 
will relate it to the above-mentioned streaking delays. 



2 Attosecond pulses of light 

It it well-known that short light pulses can be described as wave packets, i.e. 
as coherent superpositions of monochromatic light waves [SS] . The propagation 
of such wave packets can then be understood in terms of the phase-shifts of 
the individual monochromatic waves, which will be discussed in Sec. |2.1| An 
attosecond XUV pulse can be strongly stretched in time by the propagation 
through a dispersive medium, but conversely, it is also possible to compress a 
stretched pulse using a medium with negative (anomalous) dispersion. A useful 
theoretical tool adapted to describe these phenomena is the so-called Stationary 



Phase Approximation (SPA), which is introduced in Sec. 2.2 and applied to the 
propagation of attosecond pulses. As we will show later, it applies equally well 
to the propagation of photoelectron wave packets. 

We start our analysis of attosecond XUV pulses from the Maxwell equa- 
tions in a dielectric medium [TSl [S5] , and recast the wave equation into a time- 
independent form by exploiting the orthogonality of the Fourier integral expan- 
sions of the time-dependent fields 



V 2 E+^rE = -I^P-Iv(V-P), (7) 
where the electric field is expanded as a superposition of monochromatic waves 
E(t,r) = ±-Jdu J E(u;,r)e-^ t , (8) 
and the polarization is given by 

P(t,f f ) = ^JduP(u},r^e- iut . (9) 

where the fields in space-frequency are written without a tilde: E(uj, r). Finding 
solutions to Eq. Q requires knowledge of the polarization terms, P(oj,r). The 
polarization induced by a weak attosecond XUV pulse in an isotropic medium 
can be approximated by the linear response: 

P(LJ,r) = eox {1) (u,?)E(u,r), (10) 
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where the linear susceptibility, x (u, 1 ?), may vary in space due to changing 
density or composition of the medium. Strictly speaking, the medium may ex- 
hibit a more complicated evolution in time, e.g. due to ionization of the species 
constituting the medium or due to non-linear polarization, but we will not con- 
sider the latter effects here [T^l E>D] , due to the fact that the interaction between 
pulses with moderate intensity in the XUV range and matter are typically per- 
turbative. We mention, however, that such contribution may become important 
when studying the propagation of intense XUV or X-ray radiation pulses, such 
as the ones generated by XFEL facilities, with peak intensities well above the 
atomic unit of intensity, i.e. in the range beyond 10 16 W/cm 2 . 



2.1 Phase propagation of light in the linear regime 

In order to describe the propagation of linearly polarized, attosecond light pulses 
in a medium, we insert the following ansatz: E(w, r) oc exp [i<j)(oj, x)] z, into 
Eq. Q to obtain a differential equation for the phase of monochromatic waves: 



dx 2 \dx J c 2 c 2 

For slow variations of the optical properties in material, we assume \d 2 <fi/dx 2 \ <C 
\d<fi/dx\ 2 , which leads to a simpler differential equation: 

The phase of the electric field can be directly integrated as 

/x i pX 

dx' ^l + x {1) {x') = - / dx' n{x'), (13) 
-oo * ^ J — oo 

where n(x) is the local refractive index of the medium. The asymptotic phase- 
shift obtained after propagation through the material can be written as 

8{uj) = Urn - / dx' [n(x') - 1] , (14) 
x^oo c J_ oa 

which will have a finite value assuming a finite material surrounded by free 
space, e.g. that n{x) = 1 for large enough \x\. Thus, a frequency-independent 
refractive index leads to a linear phase-shift in frequency. This linear phase-shift 
causes a simple time-delay, r = S/uj, of the attosecond pulse in Eq. Q, 

Ef ree {t~T,x) = ^- [ duj Ef ree (oj, x)e~ wt exp[iwr] , (15) 

-> * ' ^ J v ' " v ' 

Shifted wave packet Free propagation Phase— shift 

compared to free-space propagation, corresponding to the undistorted waves 
Ef ree (uj,x). Eq. (151 is a direct application of the well-known shift-theorem 
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of Fourier transforms. Furthermore, if the medium is absorbing, the asymp- 
totic phase will be a complex quantity, with its real part leading to the actual 
phase-shift and its imaginary part leading to an exponential reduction of the 
amplitude, exp[i<5] = exp[i3t(<$)] exp[— !J(<$)]. 

In order to illustrate the shift of the wave packet, we compute the asymptotic 
phase acquired by an attosecond pulse passing through a material of length, L, 
with a constant refractive index, 

6 = ^[n-l]L= %[n-l]L « f X «L, (16) 
c A 2c 

where A = 2-kc/uj is the wavelength of the light wave in vacuum, and where we 
assumed that the linear susceptibility is small. The corresponding delay of the 
pulse is 

r =J^«*^, (17) 
c 2c 

where we find an intuitive linear scaling with both the length on the material 
and with the change in refractive index. |59j . 



2.1.1 Dispersion and group delays 

So far we have assumed that the linear response is independent of the frequency 
of the light. This is an adequate approximation for narrow-bandwidth laser 
pulses, but it is not a good approximation for attosecond pulses that have a 
large spectral bandwidth [BO]. In a more realistic model, the material will be 
dispersive, i.e. the refractive index will be frequency-dependent, n{oj,x) [6111621 
[63] . After passing such a material, the asymptotic (spectral) phase, 5(oj), may 
exhibit a non-linear frequency dependence, which implies that the wave packet 
in the time domain, 

E(t,x) = ^ f duo E free (uj , x)e- 1 ^ exp[i6(w)], (18) 

Free propagation Phase— function 

may be shifted and deformed with respect to the free propagation. It is conve- 
nient to expand the spectral phase in a Taylor series [59] . 

(w-w ) B , (19) 

around the central frequency of the pulse, ujq. In many practical cases it is 
sufficient to consider the first few terms in this expansion. Note that we have 
already seen that the zero-order term determines the overall phase of the pulse, 
and that the first order term determines the delay of the pulse. In the optical 
regime, most materials have a positive (normal) dispersion, which means that 
the refractive index increases as a function of frequency. This implies longer 
delays for higher frequency pulses passing through the same medium. Close to 
resonances or close to the ionization threshold of a material, the refractive index 
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may have a negative (anomalous) dispersion. The Group Delay (GD) is defined 

as 

85 

tgd(lo) = ^, (20) 



and it represents a direct generalization of the shift theorem in Eq. ( 15 ), stating 
that the derivative of the spectral phase at any frequency u corresponds to a 
delay in the propagation of the components of the wave packet in that spectral 
region. The GD is a rather abstract concept, but a useful interpretation is that 
it can be related to the time when the coherent superposition of all wave packets 
in the neighbourhood of uj is constructive. The GD is not equal to the time when 
the instantaneous temporal frequency, &rg[E(t)], equals w for a wave packet 
in general. However, we will show in Sec. |2.2| that the GD does equal the time 
of this instantaneous frequency for wave packets that are stretched in the time 
domain. 

In order to illustrate the GD concept, we now consider an XUV wave packet 
after travelling through a plasma of length L, where 



is the refractive index of the plasma, lo p — y^Ne 2 /(com), is the plasma fre- 
quency, and TV the concentration of free electrons in the plasma )64| . Using 



Eq. (161, the asymptotic phase becomes negative 



<5» = - c [n p {u)-l]L « (22) 

which implies that the phase- velocity, v p \- iase = c/n{w,x), of the monochromatic 
field is faster than c (superluminal) , because the local refractive index is smaller 
than one! The corresponding GD is, however, positive 

dS p lWp L 
tgd = ~ ~ o> (23) 

which shows that any light pulse is delayed by the plasma, and that it travels at a 
speed less than c. These effects prove to be practical for phase-matching of non- 
linear processes, where control of the superluminal phase velocity is achieved by 
simply increasing or decreasing the intensity of the laser field. This is because 
it indirectly alters the density of the plasma due to a changing ionization rate, 
which affects in turn the electron density in the plasma. 



2.1.2 Attosecond pulse trains 

Our discussion of delays is very general and it clearly applies to the propaga- 
tion of both SAP and APT. APT arise naturally from the HHG process when 
the driving laser field is many periods long. The train contains two identical 
attosecond pulses per period of the fundamental field, with electric fields hav- 
ing opposite signs. In this context, uncompensated dispersion will lead to a 
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temporal overlap and interference between consecutive attosecond pulses in the 
APT, when the individual pulse duration increases beyond a half-period of the 
fundamental laser field. An ideal APT can be written as a Fourier sum, i.e. as 
a discrete version of Eq. (fsl) , over all high-order harmonics 



E(t) = \ E 2g+i\ exp [iS 2q +i - i(2q + l)ut] 



(24) 



where 2q + 1 labels the odd harmonics that are generated from HHG. We write 
"ideal" to indicate that the APT in Eq. ( 24 ) corresponds to an infinite number 



of pulses in the time-domain, while an actual APT will have a finite duration 
comparable to that of the fundamental IR laser pulse [65, 66 . A short APT is 
likely to exhibit strong pulse-to-pulse variations. The subcycle behaviour of the 
ideal APT can, however, be understood as an average subcycle pulse structure 
in a long APT composed of many attosecond pulses. 



2.1.3 Compression of attosecond pulses using metallic foils 

As a general rule, attosecond pulses are generated with a positive intrinsic chirp 
from the HHG process [23]. However, they can be compressed in time using 
thin metallic foils in the setup depicted in Fig. IT] (a) [57] [551 E2] ■ An example of 
such XUV pulse compression is shown in Fig. [6] where the experimental data 
corresponds to propagation through a 400 nm thin Al-foil. 



p 1 
I 05 



1 



ia 



Compressed: 160 as 
Uncompressed: 180 as . 




2000 2500 



(b) 








l/^v^^x Uncompicssed^V 1 1 





500 1000 1500 2000 



2500 









^^^^^^bmpressed "^^^ 



eg 



500 



1000 1500 2000 

Time (as) 



2500 



17 



Figure 6: Attosecond pulse compression: (a) Attosecond pulses generated in Ar 
before (grey curve) and after (black curve) passing through a 400 nm Al-foil. In 
addition to a delay of the pulse of ~ 250 as, the pulse shape changes due to the 
foil, (b) The uncompressed pulse is 180 as and it is asymmetric with small post 
pulses, (c) After the foil, the compressed pulse is 160 as and more symmetric. 
The CEP has been set to be "cos-type" and the IR laser field (dashed curve) 
is plotted for comparison. Amplitude effects due to the foil are neglected for 
simplicity. 

The Al-foil provides negative dispersion in the low-energy part of the trans- 
mission window, while the high-energy part provides positive dispersion. In ex- 
periments, this property of Al-foils is used to generate close to Fourier-limited 
attosecond pulses by compensating for the intrinsic positive chirp of the HHG 
process by the negative dispersion provided by the low-energy part of the metal- 
lic transmission window. Notice the structure of the APT, with two pulses per 
period of the fundamental laser field. We mention again that the electric field 
of the two adjacent pulses have opposite phases, i.e. that they are 7r-shifted. 



2.2 Stationary Phase Approximation 



The SPA is an essential theoretical tool in attophysics because it can be used 
to evaluate integrals with complex valued integrands, such as the ones encoun- 
tered when computing time-frequency Fourier transforms. It has proven itself 
especially successful for understanding the HHG process [231 HI] , and it remains 
a useful tool for gaining better understanding of strong field dynamics, see for 
instance Ref. [701 [7H [72] . In the next subsections: Sec. |2.2.3| and|2.2.4 
apply the SPA to evaluate XUV wave packets. Further on, in Sec 



4.4 



we will 
we will 

apply the SPA to the propagation of electron wave packets ionized by attosecond 
pulses. 

In simple terms, the SPA is related to the well-known integral of a Gaussian 

(25) 



/oo 
duj exp [— /3w 2 ] 



which has a finite value provided that (3 > 0. As an example, we plot the 
integrand and integral for /3 = 1 in Fig. [7](a) and (b) respectively. The integral 
value in (b) approaches y/n rj 1.7725 as expected from Eq. (25 1. 



18 




Figure 7: Gaussian integrands and integrals: (a) A real Gaussian function 
(f3 = 1). (b) The real Gaussian integral equals y/n ~ 1.7725. (c) A complex 
Gaussian function (f3 = —i). Notice the fast complex oscillations for large values 
of |w|. (d) The complex Gaussian integral equals exp[i7r /4]. 



2.2.1 Method of steepest descent 

In more detail, the SPA belongs to a class of techniques of wide use to compute 
integrals with integrands containing the exponential of a complex analytical 
function, f(z) = u(x,y) + iv(x,y), with z = x + iy. It is closely related to 
the so-called steepest descent method which deals with integrals of the following 
form [73]: 

(26) 



K = <j) dz exp {f(z)] 



where the contour C is chosen so that the ends of the path do not contribute 
significantly. In the vicinity of a complex saddle-point z , where f'{z ) = 0, the 
complex function can be expanded as 



1 



f(z)^f(z Q ) + -(z-z ) 2 f(z ). 



(27) 



Inserting the expanded function into the exponential yields the approximate 
integral 



K 



dz exp 



c 



1 



f(z ) + -(z-z f f"(z ) 



(28) 
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where the zeroth-order, f(z ), and the second-order, f"(z ), coefficients are 
complex. In order to evaluate this integral, we parametrize the contour: z(p) = 
zq + pe lip , with p being a real parameter and ip being the angle at which the 
contour crosses the saddle-point at z in the complex plane. After this change 
of variable we get 



K « exp [/Oo) + ip] 



dp exp 



^ exp[z2^]/"(z ) 



(29) 



where the integral runs over the real parameter, p, but the integrand remains 



complex. Using Eq. ( 29 ), it is, however, simple to choose an angle for the contour 
so that the exponent becomes negative and real, 2ip + arg[/"(zo)] = tTj which 
turns the integrand into a simple real Gaussian with j3 = |/"(zo)|/2 in analogy 



with Eq. (25). A technical point is that we must have — n/2 < tp < n/2 for the 
contour to run in the "right" direction, i.e. with the real part of z increasing, 
but this is typically the case in physical problems. The final result simplifies to 
the following beautiful expression 

K a iex P [/O )]^^, (30) 

which we can interpret as the "height" exp[/(zo)] times the "width" ^2n/ f" (z ) 
of the integrated function close to the saddle-point. 



It is now of interest to compare the simple integration of a Gaussian, Eq. ( 25 ) 



with the result of our complex analysis, Eq. (pOft, and conclude that the in- 



tegrand in Eq. ( 25 1 may be complex. The integral converges provided that 
5R(/3) > 0, but also in the special case of Sft(/3) = 0, provided that 9(/3) ^ 0. 
In the latter case, the convergence can be understood as being due to rapid 
complex oscillations for large values of |w| as shown in Fig. [7](c), leading to a 
finite complex integrated value in Fig. [7](d) of y/w exp[wr /4]. 

2.2.2 General equation for SPA 

Inspired by the results above, we now consider a typical frequency domain in- 
tegral of a more general function: E(u>) exp(if(oj)), where E(u) is a smooth 
amplitude function and f(u) is a smooth phase function. The phase function 
may contain several saddle-points, also called "stationary" points. In regions 
where f(ui) is varying fast, we expect the contribution to the integral to be 
small in analogy with Fig. [71(c) and (d). Conversely, we expect significant con- 
tributions from regions close to the stationary frequency, Ui s , where the first 
derivative of the phase function vanishes: df /doj\ OJs = f'(u) s ) = f' s = 0. Then, 



using the truncated Taylor expansion in Eq. (271 in the vicinity of the station- 



ary points and replacing the exact integrand by the approximate local Gaussian 
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function, one gets the final result for the SPA approximation [75] : 



/ / 2tt{ 
du> E(u)exp[if(u)] « J^.J — exp [*/.], (31) 
-oo D \ Js 



where / s = /(w s ), /" = d 2 //<9aj 2 |^ and B s = E(uj s ). The sum runs over all 
the roots of the stationary-phase equation: 

t; " °' (32) 

where the roots w s are labelled by the index s. The number of solutions depends 
on the properties of the phase factor, f(ui). In the simplest cases, which we will 
consider here, there is only one solution. The SPA is derived assuming that the 
second-order term is dominant as compared to higher-order phase terms, e.g. 
that (oj— w s ) 2 /"/2 3> (u>— w s ) 3 /"'/6, in the neighbourhood of each saddle point, 
\w — lj s \ < 7T. If the second-order phase term happens to be zero, f"/2 = /3 = 0, 



this would correspond to an artificial divergence in Eq. (31 ). This divergence is 
most likely non-physical [and it may happen close to the cut-off in the HHG] and 
the result would be finite if the third-order phase term in the Taylor expansion 



of the exponent was used in the integral. Having derived the SPA, Eq. (31 1, we 
now need to demonstrate that it can help us to make physical arguments about 
the propagation of XUV wave packets. 

2.2.3 Stretched XUV wave packets 

It is well-known from ultra-fast optics that uncompensated dispersion can lead 
to stretched pulses in the time domain [8]. It is illustrative to use the SPA to 
calculate the properties of such stretched pulses. Consider an attosecond pulse, 



as expressed in Eq. (18), at x = 0, with a rather broad Gaussian envelope, 

E(u) = E exp[-g{u> - lu ) 2 }, (33) 
and with a phase factor f(ui) = —ut + S(lo), where the spectral phase is 

5(u) k6 + (uj- u> )5' + i(w - wo)X. ( 34 ) 



with a large quadratic chirp, |<5q|/2 ^> g. The SPA expressed in Eq. (311 then 
implies that the magnitude of the amplitude factor will be dominated by the 
contribution of a local value, E s = Eq exp[— g(u> s — w o) 2 L occurring at the 



stationary frequency, ui s , as determined by solving Eq. (32 ). This is quite similar 
to the way in which a delta function operates, selecting a specific value of a 
function. Here it arises from the small constructive window in the complex 
integral that occurs around the stationary frequency. 

The next task is to find the stationary frequency, uj s , which can be seen as 
the frequency that makes the dominant contribution to the wave packet in a 



given point in space-time. The stationary-phase equation, Eq. (32 1, becomes 
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first-order: f'(uj) s» S' + (co — u>o)^o —t = 0, which we use to find the stationary 
frequency for a given time, t: 



1 , 

w o + (t 



6' n % 



(35) 



We can rewrite Eq. (35 ) to see that, due to the chirp induced from propagation, 
any frequency of the pulse, lj = ui s (t), "arrives" at a unique time, t = 5' + 
[to s — (JqWo ~ 3' s > wm ch equals the GD as defined in Eq. (20 1. Furthermore, 



one can show that the stationary frequency equals the instantaneous temporal 
frequency of the wave packet: 



U) a (t) 



d_ 

Ft 



arg {£(*)} 



(36) 



which provides a simple link between time and frequency for stretched pulses. 
In the derivation above, it is only required that the spectral envelope is "broad" 
as compared to the spectral width induced by the quadratic chirp. Any other 
broad envelope, E(u), will also work, leading to a mapping of the intensity from 
the spectral domain to the temporal domain, 



Temporal 



1 

M 
Constant 



Spectral 



(37) 



because the stationary frequency increases linearly in time as shown in Eq. ( 35 ) 



The mapping of the spectral envelope to the temporal envelope is scaled by the 
inverse spectral chirp, 1/8q, so that an increased spectral chirp leads to a more 
stretched pulse in the time domain. At the end of Sec. |4.4[ we will see that 
the properties of these stretched light pulses are similar to the way in which 
non-relativistic electrons are dispersed through propagation in vacuum. 



2.2.4 Ultra-short Gaussian pulse 



Consider again an attosecond pulse, Eq. (31 ), with a Gaussian envelope, Eq. (33 1, 
but this time the envelope is not "broad" compared to the quadratic spectral 
chirp, instead, as in the cases of interest in attoscience, one have XUV pulses 
so that g \6q |/2. Then, for the SPA to be meaningful, the Gaussian envelope 
must be included in the phase factor so that it becomes complex: 



f(u>) = S(u>) — ojt + ig(uj — ujq) 2 



(38) 



Using the asymptotic phase to second order, Eq. (34 1, the stationary-phase 



equation, Eq. (32), remains first-order, with complex and time-dependent coef- 
ficients 

/'(«) *6' + SH(u - wo) - * + *2 g (w - lu ) = 0, (39) 

with the root: 

i2g 



W a (t) =UJ + JJJ 1 - 



(t-S' ). 



(40) 



22 



We mention that it is common to obtain complex stationary points in the SPA 
in the case of steep variations of the spectral amplitude. 

The real part of w s is identified as the instantaneous frequency, as was done 
in Sec. 2.2.3 while the imaginary part of lo s limits the intensity duration (or 
probability distribution) of the attosecond pulse in time. At the origin, x = 0, 
the intensity of the attosecond pulse will be large only at a time, t, close to 
the group delay: t m S' 0l where the imaginary part of the stationary frequency 
vanishes. Inserting uj s (t) into the right-hand-side of Eq. (|3T|) yields the result: 



^ )= 2^V/3 eXP 



i(5 - oj t) 



4/3 



(41) 



where p = g — iS'^/2. In this particular case, the result of the SPA is exact 
because it was applied to a Gaussian pulse, but the method is more general and 
it can be applied also to other spectral envelopes. This is achieved by replacing 
the third term in Eq. (38) by the more general expression: \n[E{uj)]/i. 



In the special case of a Fourier-limited pulse, <5q = 0, the real part of the 
stationary frequency becomes the central frequency of the pulse at all times; 
while the imaginary term remains time-dependent: Ui s (t) = ujq + i(5' — t)/2g. 
This is closely related to the Fourier-limited pulse shape, which has a constant 
carrier frequency equal to the central frequency. The imaginary part of the 
stationary frequency imposes the finite duration of the pulse in time. 



Summary of attosecond pulse propagation 

In this section, we have reviewed the basic properties of light pulses that prop- 
agate through a medium in one dimension. Attosecond XUV light pulses are 
a special kind of light pulse with a large spectral bandwidth. The larger the 
bandwidth, the more sensitive the pulses are to the dispersion induced by the 
medium. Using an appropriate medium, it is possible to correct for spectral 
chirps and reduce the attosecond pulses to their shortest possible duration. 
Once the attosecond pulses are adjusted, they can be used to ionize atoms at 
extremely well-defined times |74j . We now review the properties the properties 
of non-relativistic photoelectron wave packets. 



3 Photoelectron wave packets 

Attosecond XUV pulses have been used to initiate photoelectron wave packets 
at well-defined times, but how does such a wave packet evolve after the ioniza- 
tion event? In this section, we perform an analysis of the evolution of these wave 
packets and we shed light on the similarities with the propagation of attosec- 
ond XUV pulses, Sec. [2] Although our final goal is to describe photoionization 
as a time-dependent process, we start our study with a theoretical treatment 
of positive energy photoelectron wave packets in the presence of a static bind- 
ing potential. More details on the actual attosecond-XUV transition are given 
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in Sec. [4j while the important features related to laser-probing in long-range 
Coulomb potentials are given in Sec. [5] 

In the absence of an external field, the time-evolution of a non-relativistic 
photoelectron is described by the time-dependent Schrddinger equation (TDSE) 

m 

H |*(t)) = [T + V )\*(t)) = in^\*(t)), (42) 

where Hq is a time-independent Hamiltonian, and ^ = ^(i, r) = (r|\l/(i)) is 
the time-dependent wavefunction. The Hamiltonian consists of two parts: T = 
—h 2 \7 2 /(2m) being the kinetic energy operator, and Vq = Vo(r) being a time- 
independent, single-electron potential describing the binding of the electron to 
the atom. For simplicity, we limit ourselves to systems that can be described 
by a single electron, namely the SAE approximation. A common approach to 



solving the TDSE, Eq. (42 ), is to first find the solutions to the time-independent 
Schrddinger equation, 

HMr) = e a ^ a (r), (43) 

where e a is the energy (eigenvalue) of the solution (eigenstate), ip a (r) labelled by 
a set of quantum numbers a. In the following, these time-independent solutions 
are referred to as states. They form the basis that spans the complete space in 
which the electron moves. 

For a spherical potential, V(r), the states can be written on a spherical 
basis, ipa(f) — Rn,\{r)Y\, m {9,(j)), where R n ,\ is the radial wavefunction and 
Y\,m is the spherical harmonic, for quantum numbers: a — [n, A, m], with spin 
neglected for simplicity. It is possible to reduce the Schrddinger equation into 
separate radial equations, each corresponding to a specific angular momentum 
quantum number, A. To emphasize the similarities and differences between the 
one-dimensional propagation of light wave packets [Sec. [2], the effective radial 
electron states, u n .\(r) = rRn,x{r), are used. These effective radial states can 
be interpreted as the solutions of a one dimensional Schrodinger equation for 
r > with an effective potential 

[T + V r (x) + V x {x)}u n , x (x) = e a u n . x (x), (44) 

V (x) 

where T = — h 2 /(2m)d xx is the radial kinetic energy, V r (x) = —e 2 /(4ireox) is 
the radial potential (explicit for hydrogen) and V\(x) = A(A + l)h 2 /(2mx 2 ) 
is the centrifugal potential that depends the angular momentum through A. 



Starting with Eq. (44 ), the radial variable r is relabelled x to ease the comparison 
between the propagation of one-dimensional light pulses and of photoelectron 
wave packets. The boundary condition for all states is u n ^\(0) = 0, and the 
states are taken to be real A photoelectron wave packet, in angular subspace A, 
can then be can be written as an integral-superposition of continuum states, 

/>oo 

^\(t,x)= / de a(e,t)u e! \(x) exp [-iet/h] , (45) 
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where a(e, t) is the complex amplitude for the energy-normalized radial state, 
u €t \(x), at energy e > . 



3.1 Phase propagation of photoelectrons 



The propagation of a photoelectron wave packet, Eq. (45), is governed by the 
way in which the phases of the continuum states u t ^\(x) vary as a function of 
energy, e. The Wentzel-Kramers-Brillouin (WKB) approximation is well suited 
for studying continuum states within a semi-classical framework [75] . The WKB 
approximation relies on a similar ansatz as in Sec. |2.1| for XUV light. For an 
effective potential which vanishes at infinity, lim^^oo Vq (x) = 0, the effective 
radial wavefunction is taken to be 

u e (x) cx 9 {exp[i0(e, x)]} — sin[0(e, x)], (46) 

where reference to A is dropped for compactness. Inserting the WKB ansatz 



into the Schrodinger equation, Eq. (44), leads to a differential equation for the 
phase of the state, 



h 2 d 2 <\> h 2 

2m dx 2 2m V dx 



e-V (x). (47) 



Further, assuming that the local momentum of the photoelectron varies slowly, 
\p(x)/h\ = \d(j)/dx\ ^> \d 2 <j)/dx 2 \ 1 / 2 , the above equation can be simplified, 



£ = -y/2m[e-V (x)}, (48) 



and then integrated, 



(e,x,x ) = - I dx' y / 2m[e- V (x')] 



dx' p(x'), (49) 



•I'll 



where p(x) is the local momentum of the electron, and xq is from where the 
electron "starts" . In applications, the parameter xq is associated to a classical 
turning point, see below. The local momentum corresponds to the local kinetic 
energy, which is the total energy minus the local potential energy. The phase 



of the wavefunction varies faster as the local momentum is increased. Eq. ( 13 ) 



and ( 49 ) are now compared to highlight the differences between the propagation 
of XUV light and photoelectrons. In a medium with negligible dispersion, the 
phase of light is approximately proportional to the photon momentum, w/c, 
leading to a simple translation of the wave packet in space-time according to 



the shift theorem, Eq. ( 15 ). Similarly, in regions of space where the action of the 
potential is negligible, Vq(x) -C e, the phase of the electron is proportional to 
the electron momentum, but this corresponds to the square root of the energy 
and frequency, p(x) cx y/e, in contrast to the linear frequency dependence of 
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light. This implies that the electron wave packet will not maintain its shape 
under propagation. This non-linear phase dependence leads to a broadening 
effect referred to as "quantum diffusion" and it originates physically from the 
fact that electrons of different energies travel at different speeds. In the more 



general case with dispersion, we may rewrite Eq. ( 49 1 and define an "electron- 
susceptibility" — —Vo(x) /e in analogy with the linear susceptibility phase of 
light, x^- The electron phase is real as long as e > Vq, but it will be imaginary 
if e < Vo, which implies that the electron must tunnel through some barrier and 
that the transmission amplitude will be exponentially damped. 

When writing the WKB solution, it is natural to include \d 2 4> / dx 2 \ to leading 
order [75] ■ This leads to the well-known form for the WKB state 



(WKB) 1 

u , , , cx — ; exp 



±1 [ dx'ptx') 



(50) 



where the two complex solutions (±) describe an outgoing and incoming elec- 
tron, respectively. The probability density, p(x) — \u^± B \x)\ 2 cx l/|p(x)|, 
is smaller where the electron moves rapidly. Physically, we understand this 
because the electron will spend a shorter time in a place when moving with 
a greater velocity. Alternatively, we can form an energy normalized and real 
WKB solution, 



S e (x,x ) + 4> Xo 



(51) 



which is valid for x > xq, where x — xq is a classical turning point, i.e. it 
separates a classically forbidden region, e < Vq(x) for x < xq, and a classically 
allowed region, e > Vq(x) for x > xq. The phase of the state is a sum of the 
semi-classical action, 



\s t (x,x ) = \ f dx'p(x'), (52) 

I X 



plus the quantum mechanical reflection phase, <p Xo . Depending on the type of 
potential and on the energy of the electron, the quantum mechanical reflection 
at the classical turning point can be regarded as either "hard" or "soft" [75]. 
Photoelectrons in a pure Coulomb potential, with no angular momentum, will 
suffer a hard reflection at xq = 0, with a reflection phase of <f> Xo — 0, so that the 
wavefunction is exactly zero at the reflection point, i.e. u^ KB \o) cx sin(0) = 
0. Photoelectrons with angular momentum, A > 0, will reflect in an effective 
potential composed of the Coulomb potential and the centrifugal potential at 
Vo(a;o) = V r (xo)+V\(xo) = e. This combined barrier typically leads to reflection 
points further out from the nucleus, xq > 0, and to softer reflections. A soft 
reflection implies that the photoelectron wavefunction can extend a bit into the 
classically forbidden region, x < xq. The soft reflection limit, i.e. when the 
effective potential is slowly varying compared to the wavelength of the electron, 
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can be calculated using a linear potential so that the wavefunction becomes 
an Airy function with a reflection phase of 7r/4 [75], but this is typically not 
a good approximation for photoelectrons that we consider here. An accurate 
determination of the reflection phase, <f> Xo , requires a fully quantum mechanical 
theory beyond the WKB approximation and ab initio calculations of such phases 
present a challenge to theory for large atomic systems or complex molecular 
systems. 

3.2 Quantum diffusion and Wigner delays 

As an example, consider the situation illustrated in Fig. [HJ where an electron 
passes through a short-range attractive potential. 
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Figure 8: Electrons passing through an attractive potential, Vq(x). The real 
WKB solutions (thin black curves), Eq. (51), are compared to the real free 
electron states (thick grey curve) at three different energies: e = 1,2,3 atomic 
units of energy (27.2 eV). The phase of the electron wavefunction varies more 
rapidly in the potential, which leads to an accumulated phase difference, <5, 
compared to the free case. (Note that the modulations of the wavefunctions 
should be interpreted in the third dimension of the graph, and not as an energy 
modulation.) 

The electron is classically allowed to pass through the potential, but as it does 
so, it will acquire a quantum phase. The asymptotic phase difference, between 
an electron wave propagating through a short-range potential, Eq. (49), and 



free wave propagation, is defined in analogy with Eq. ( 14 ) 



6(e) 



1 



lim 



dx' \p(x) - p ] 



(53) 



where po = \j2vne is the free (asymptotic) momentum of the electron. The GD 
concept, Eq. pOj), can also be applied to electron wave packets using Eq. (53): 



. > .9(5 as 



(54) 
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The electron GD then results only from the interaction with the potential, since 
the intrinsic delay due to quantum diffusion has been subtracted. Using the 
asymptotic phase to calculate the delay of electron wave packets was first done 
by Eisenbud, Wigner and Smith [77J EH ES] , hence the name Eisenbud-Wigner- 
Smith delay, (or Wigner delay for short). In the simple scattering example given 
above, the asymptotic phase shift equals twice the radial phase shift: once for 
the incoming electron (0 > x) and once for the outgoing electron (0 < x). 
Regarding photoionization, it can be considered as a half- collision because the 
electron only moves out of the atom (0 < x). As compared to free propagation, 
the delay of the wave packet in photoemission is, therefore, the derivative of the 
radial asymptotic phase-shift. 

For a weak potential, \ Vo\ <?C e, the local momentum is p{x) ~ p [1 — V~o(x)/2e], 
which leads to an asymptotic phase 

„, < 1 frn f 00 , , 1 pm 

s U*-ni2eL dxVoix ^-nV^ Iv > (55) 

where the potential integral can be written as Iy ~ eyLy, in terms of an 
effective "height" ey and "range" Ly of the potential well. As expected, in 
this lowest-order approach, the phase acquired by the electron increases linearly 
with both the height and range of the barrier. The Wigner delay is 

l™ Iv ~ fm ty 
T * = V8^~V 8^ iv ' (56) 

It takes a shorter time, t\ < 0, for an electron to pass a weak attractive potential, 
ey < 0, than through free space. Physically, this is due to the higher local 
velocity of the electron in the potential valley. The opposite is true for a weak 
repulsive potential, where the electron is slowed down on the potential hill. 
Furthermore, the timing of slow electrons is affected more than the timing of 
fast electrons, due to the longer total time that the slow electrons spend in the 
potential. In terms of numbers, this simple analysis leads to a delay of 5.9 as per 
eV and nm of the potential barrier, given an electron with one unit of atomic 
energy (27.2 eV). 



3.3 Phases and delays in Coulomb potentials 



Photoelectrons created from neutral atoms are submitted to the long-range 
Coulomb potential of the remaining ion, Vc = —C/x, where the constant is 
C = Ze 2 / (Aireo) with Z = 1 for a singly charged hydrogen ion. In order to 
study this long-range interaction in detail, we consider the WKB approxima- 
tion, Eq. (15 1 1) . Then, the local momentum can be expanded as 
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with po = \j2mt being the asymptotic momentum in the remote region of space 
where e > \C/x\. Using Eq. and ([57]), the total phase of the real WKB 
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state is asymptotically 

S e (x,x ) 



kx + — \n{x) — kxo — — ln(xo) 



= kx 



\n(2kx) 
kao 



T 



(58) 



where k — po/n is the wave number; ao is the Bohr radius; and Z = 1. Line 
1 of Eq. ( 58 1 corresponds to the WKB notation with an approximate reference 



to the reflection point xq, while line 2 follows standard notation from atomic 
physics with reference to the absolute scattering phase, r\ = rjk,\- The real WKB 
state varies as 
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(59) 



for x — > 00. Interestingly, the phase does not settle into free particle behaviour, 
instead, the phase diverges logarithmically! What does this space-divergent 
phase mean for the delay of the photoelectron wave packet? When applying the 
frequency derivative to the the logarithmic term, we find a delay that depends 
on the spatial position of the wave packet, 



T~Ln(k,x) = h 



We 



ln(2fex) 
kao 

nook* l M^T T 



(60) 



where we have replaced this spatial-position of the electron, x, by the approx- 
imate classical position, x c i(t) s=s vt = hkr/m, with r = t — to being the time 
that has passed since the electron was ionized by the attosecond pulse at t = to . 
Interestingly, also diverges logarithmically as r increases and it cannot be 
neglected. As a consequence, Tm will completely dominate over any (short- 
range) Wigner delay, t\ -C tlr <C t, as r — > 00. This implies that there is 
no absolute delay of the Coulomb wave packet with respect to free propagation. 
To define a physical delay between two photoelectron wave packets of different 
energy is already difficult for free propagation, due to the different velocities of 
the photoelectrons as discussed with Eq. (54 1, but it is even more difficult for 
Coulomb wave packets due to the spatial dependence of the phase. In the special 
case of two different Coulomb wave packets of the same energy, the logarithmic 
delays cancel and their relative delay then equals the difference in Wigner delay. 
In this way, we can talk about the delay of a Coulomb wave packet, but only 
relative to a Coulomb reference, e.g. the hydrogenic system. To summarize 
this analysis, for the propagation of photoelectron wave packets, we find that 
the approximate position is predominately given by the linear, free propagation 
relation, pox — er — 0. The logarithmic, long-range phase is of intermediate im- 
portance, and it dominates over the absolute, short-range phase. On the other 
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hand, only the short-range phase contain detailed information about the atomic 
potential. 

As already mentioned, the absolute determination of the short-range asymp- 
totic phase-shift, 77, is a difficult problem, because it depends on the detailed 
atomic potential, which is modified by electron-electron interactions at short 
range. In the special case of hydrogen system of nuclear charge Z, and for a 
given wave number, k, and angular momentum, A, the phase is known analyti- 
cally [75], 

ff fc| A = arg{r[A + l-tZ/(fcao)]}, (61) 



(H) 



with r<» = / °° dt t z ~ x exp[— t] being the complex Gamma function. This result 
can be generalized for any spherical atom, 



Vk,\ = 0"fe,A + 6k, A, 



(62) 



where 6k,\ is the phase difference compared to hydrogen that is induced by the 
short-range deviation from the pure Coulomb potential. The WKB states are 
compared with the exact Coulomb states for hydrogen in Fig. [9] 
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Figure 9: Exact hydrogen wavefunctions (thin black curve) for A = 0, I and 2, 
corresponding to s-, p- and d- waves in the continuum at three different energies 
5, 15 and 25 eV. The asymptotic approximation (grey thick curve) is shown for 
comparison. It will be explained further in Sec. [4] and [5] but we mention here 
that it utilizes the long-range WKB phase matched to the correct asymptotic 



phase-shift, while the amplitude prefactor is taken to be constant [see Eq. ( 46 )] . 
(Note that the modulations of the wavefunctions should be interpreted in the 
third dimension of the graph, and not as an energy modulation.) 
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The exact analytical states of the hydrogen atom are called the regular Coulomb 
functions of the first kind [50] EI] • The WKB states converge quickly to the 
exact solution, if the correct scattering phase is known and applied [76] . 



4 Photoionization on the attosecond timescale 

In this section, we will discuss, within the framework of lowest-order time- 
dependent perturbation theory [52] [S3], the photoionization of atoms in the 
presence of attosecond XUV pulses (SAP or APT). Photoionization is a funda- 
mental process in light-matter interaction, which has been studied extensively 
using time-independent methods, in parallel with the implementation of syn- 
chrotron radiation sources, see for instance the seminal papers about Cooper 
minima and photoelectron phase-shifts [5jJ [55] . Here we start by studying pho- 
toionization in the time-domain, but as we proceed we shall recover the more 
traditional, spectral representation for the determination of the phase-shifts of 
the photoelectron states, from which one can recover relative delays between 
photoelectron wave packets. 

Although attosecond pulses of high-order harmonics have relatively high 
peak intensities, as compared to conventional sources of XUV radiation, single- 
photon absorption remains the dominant mechanism for atomic ionization. Then, 
due to this fluence of XUV or X-ray photons being sufficiently small, the interac- 
tion between attosecond pulses and atoms can be approximated using first-order, 
time-dependent perturbation theory [55] [S3]- The field-free TDSE, as written 



in Eq. ( 42 ) , is satisfied by the time-dependent field-free states of the atom, 

|a(t)) = | a) = | a) exp[— iu) a t], (63) 

where a labels the quantum numbers for a given energy, e a — huj a . The inter- 
action between the XUV field, Vj(z,t), and the unperturbed, initial state of the 
atom, \i(t)) with energy = TiUi and a binding energy of I p = \ei\, leads to the 
creation of a photoelectron wave packet 

!*(*)> « |?(*)> + Jt a ^(t) |a(t)>, (64) 



where the first-order complex amplitudes are |82j 

a«(<) = i f dt' <a(0|VKM')RO>- (65) 

In the case of photoionization, the interaction Hamiltonian, Vj(z,t) = ez E(t), 
is the dipole-interaction operator for the attosecond XUV field, here given in 
length gauge for linear polarization along z. For the case of short, coherent 
pulses, e.g. attosecond pulses, the interaction can be expanded as 

Vi(z,t) = ez -!- / dCl' E(n')exp[-iQ% (66) 
2tt J 
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where B(fi') are the electric-field, spectral components of the attosecond pulse. 
After the interaction is over, when the attosecond pulse has passed, the first- 
order complex amplitude, or S-matrix element (not to be confused by the elec- 
tron action), is proportional to both the dipole matrix element, z ai = (a\z\i}, 
and to the Fourier transform of the electric field, 



t— > oo 

e 



z a i J dt E{t) exp [i(u> a — tOj)i] 
w E(Sl), (67) 
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at the energy-conserving frequency, 51 = uj a — Wj, where \a) is a continuum 
state under the single-active electron (SAE) approximation. This shows that 
the Fourier components of the attosecond pulse are continuously mapped on the 
complex transition amplitude associated to the photoionization process. Most 
attosecond pulse-characterization schemes take advantage of this mapping from 
XUV light to photoelectron states, see for instance Ref. [55] and [32] [33J for the 
attosecond streak camera and the RABITT method respectively, which we will 
discuss further in Sec. [5] 

More generally, in order to determine the outcome of the interaction, e.g. 
the momentum distribution, the complete photoelectron wave packet can be 
projected on a momentum state, 

S %/i = Hm $*(*)>, (68) 

of the field-free Hamiltonian with energy, e% — h 2 k 2 /(2m). The real expec- 
tation values, \S^u\ , yield then the corresponding probability density of the 
momentum states. These momentum states are not plane waves, but they do 
approach a plane- wave like behaviour in the asymptotic limit, r — > oo. Their 
wavefunction, fj:{r) — (r\k), can be expanded on partial- waves as [83j 

n (r) = (Stt) 3 / 2 J2 ^ L e- ^ ^Yi M {k)Y L Mf)R k Ar), (69) 

L.M 

where the inserted scattering phases, i L e~" ,L ( k \ are designed to form an unified 
phase front, of all partial waves (L, M), in the forward direction, k. For electrons 
ejected from a neutral system, outgoing waves have the following generic form 



(as we will discuss further in Sec. 4.3) 



u k U L ] ( r ) « exp |i kr+ + vUk) - ^ j 



(70) 



where the asymptotic phases clearly cancel with the phase factors in Eq. ( 69 1 
for each partial wave. However, there remains a logarithmic term due to the 
long-range Coulomb potential, as derived in Sec. |3.3[ which shows that the 
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asymptotic behaviour of the momentum state is not exactly that of a simple 
plane wave. Furthermore, it is only the phase difference between the partial 
waves that is physical, so that the total phase of the state can be set arbitrarily. 

In this way, the spectral intensity distribution of attosecond XUV pulses 
can be determined using electron spectrometers, such as time-of-flight (TOF) or 
velocity-map imaging (VMI), to count the number of photoelectrons reaching 
the detector, and then correcting for the dipole matrix element, i.e. for the 
"quantum efficiency" [86l HZl [88] , as shown in Fig. [I] We stress, however, that 
in either case, information about the temporal structure of the XUV light pulse 
is not accessible because the phase of the complex amplitude is not measured. 



4.1 Snapshots of photoelectron wave packet 

By considering a simple "flat-top" XUV attosecond pulse, we can observe 
"snap-shots" of how the photoelectron wave packet builds up with time. The 
flat-top attosecond pulse is of a given duration, 77, and it has a central fre- 



quency of TlVL > I p . Fig. 10 illustrates the evolution of the complex amplitudes 
and the corresponding reconstructed photoelectron wave packets, ISf^ 1 '^)), for 
increasing interaction durations, 77. 
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Figure 10: Complex amplitudes (left) and "snapshots" of the corresponding 
photoelectron wave packets (right) after photoionization of increasing pulse du- 
ration, 77: (a,b) 120, (c,d) 240, (e,f) 480 and (g,h) 2400 as. (a,c,e,g) The prob- 
ability distribution over energy (black line) becomes narrower as 77 increases. 
(b,d,f,h) In space, this corresponds to an expanding electron wave packet. The 
size of the photoelectron wave packet is well estimated by the classical distance 
x — v c ti (black circle). This is a one-dimensional model calculation where the 
continuum states are plane waves and the dipole transition matrix clement is ap- 
proximated as constant. The carrier frequency of the light is Kl = 41 eV and the 
binding energy is I p — 13.6 eV, resulting in a central energy, e c — TtQ — I p = 27.1 
cV = 1 atomic unit, of the photoelectron. The photoelectron wave packet is 
complex; the real and imaginary parts being shown in red and blue, respectively. 

From standard time-dependent perturbation theory, we know that the total 
transition probability of a perturbative excitation increases linearly with time, 
P = J de a \a a (t)\ 2 oc 77, as long as the depletion of the initial state can be 
neglected. The photoelectron wave packet, in Fig. [TUJ behaves indeed as ex- 
pected since it extends linearly in space as a function of the interaction time, 
77. The spatial extent of the wave packet is approximately x ps v c Ti, where 
v c = ^/2e c /m is the classical velocity of the photoelectron after ionization and 
77 is the duration of the interaction. Inside this "classical" extent of the wave 
packet, the probability density is roughly constant. 

Because the electron is ionized at a "well-defined time" , namely during 
the brief interaction with the attosecond pulse, the uncertainty principle dic- 
tates that its energy content must be very broad. This can also be verified in 



Fig. 10 (g,e,c,a) for progressively shorter interactions. It may be of interest to 
know how long the photoelectron wave packet remains "localized" in time and 
space. At time r after the ionization, the electron wave packet will have spread 
in space to an extent of Ar, which can be estimated from a classical ensemble 
of velocities, 

Ar = (v max - V m in) T PS -T (71) 

mv c 

where v min , v c and v max are the minimal, central and maximal velocity, corre- 
sponding to an energy width of Ae at a central energy of e c = mv^/2. Further- 
more, the temporal extent can be estimated as 

Ar Ae 

At « — ps — r. (72 
v c 2e c V ' 

Interestingly, from this simple analysis, we find that the photoelectron wave 
packet will extend its duration linearly in r. For photoelectrons ionized by at- 
tosecond pulses, with a large bandwidth Ae ps e c , this implies that the temporal 
extent of the electron wave packet is approximately equal to its time of propaga- 
tion, At ps t. In other words, we should not consider these photoelectron wave 
packets as attosecond quantities. We will return to the propagation of photo- 



electrons in Sec. 4.4 for a more detailed analysis, but first we must consider the 



actual photoionization process in more detail. 
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4.2 Exponential turn-on of ionization 

In the previous subsection, we considered a "flat-top" XUV pulse shape, but 
envelopes with such sharp features represent a rather poor idealization of the 
often smoothly varying physical pulses. Rather than turning on the ionizing 
field instantaneously, as was done above, the electric field can be smoothly 
increased from minus infinity using a "slow" exponential turn-on [82j . Even 
though this represents an infinitely long interaction, we will show that it is 
extremely important for a deep theoretical understanding of photoionization 
phenomena also on the attosecond timescale. The electric field is written E(t) — 
iE exp[—i£lt + 0t], including an attenuation factor, /3 > 0, that reduces the field 
strength as t — >• — oo. The attenuation factor ensures that the probability for 
the emission of a photoelectron is zero at remote negative times. Inserting the 



attenuated field into Eq. ( 65 1 yields 



a«(i) = lim E -z m E ( 1 ) e i^ +0 t 

= \z ai E (w5{Au a ) - ip-r^—) e lA "°', (73) 

with Aw Q = uj a — uji — O, and where the limit corresponds to the "well-known 
theorem from complex function theory" |89j . The notation used requires some 
explanation: We write S for the Dirac delta function, and p is written to indicate 
that the discontinuity of 1/Aw a has been removed, so that when integrated the 
result is Cauchy's principal value. We mention that this limiting form can be 



recovered by multiplying the denominator of Eq. ( 73 ) by its complex conjugate 
and then identifying the limits of the delta function, j3 / (f3 2 + Aw^) — ► ir8(Auj a ); 
and the principal value, Aw Q /(/3 2 + Aw^) — > p(l/Aw a ). 

Similarly, we can model a slow turn-off using E(t) — i_Eexp[— ittt — f3t] for 
t > as t —¥ oo, giving the same result as in Eq. ( |73[ ), but with the opposite sign 
of the principal value part. Combining the slow turn-on and the slow turn-off, 
therefore, leads to symmetric cancellation of the principal value contribution. 
The resulting complex amplitude, corresponding to an interaction with a field 
of central frequency f2 for times going from minus infinity to plus infinity, shows 
that energy must be conserved, 

SjJJ - ez m E 2n8{e a -a -MT), (74) 



and that only an energy-conserving state, with e K = h k 2 / (2m) = hil + ej, will 
be populated. This result corresponds to a Fourier transform of a continuously 



oscillating electric field, Eq. (67 1. We stress that it is valid as t — > oo, which we 



interpret as the electric field being turned off. 



The interesting aspect of the exponential turn-on model, Eq. (73), is that it 



allows us to reconstruct a photoelectron wave packet at any finite time, t, while 
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the ionization is ongoing and the electric field is on, 



|* (1) (*)>=£«£ 1) (*) |«(*)>= 

Consermnj Principal value 

eE ' 



[n^2z ci \c) -ipjZ e Z *\ \ a ) exp[-iu K t] , (75) 



Temporal phase 



Time-independent part 
where the total phase of the wave packet is evolving at a common angular fre- 



quency, uj Ki corresponding to energy conservation. In Eq. (75), the time, t, is 
finite. This implies that all states can be populated, also those that do not 
conserve energy. In fact, the time-independent part of the wave packet consists 
of two terms: the first one is associated to the energy-conserving states, while 
the second term is a non-trivial principal-value superposition of non-energy- 
conserving states. The latter states are out of phase by ±7r/2 compared to the 
energy-conserving contribution depending on whether the energy denominator 
is positive or negative. This total phase displacement of it can be understood 
in analogy with a classical pendulum, which is set in motion by an external 
oscillating force either in an over-driven or under-driven mode. Here, because 
we are dealing with quantum mechanics, we need to sum over all possible "pen- 
dulums" corresponding to the continuum states above and below the resonance 
located at e K . 

More explicitly, the energy-conserving transition goes from the initial, partial- 
wave state, |z) = \m,£i,mi), to a superposition of energy- conserving, partial- 
waves states, |c) = |k, A c , /i c ), where the summation over c implies that A c = £i± 
1 and /i c = mi for linearly polarized light. Similarly, the non-energy-conserving 
states can be written explicitly as partial-wave states, |a) = A q ,/i q ), where 
a runs over both the bound and continuous states, but k 1 ^ n. 

4.3 Asymptotic form of the photoelectron wave packet 

It is not easy to directly infer the form of the first-order wave packet from 



Eq. (75), but we expect that it should be an outgoing wave. To see this, we 
first note that the energy-conserving part is real and that the principal-value 
part is imaginary, at t = 0. Next, we consider the asymptotic form of the radial 
wavefunction for a photoelectron state with e = e a — e K i, corresponding to 
quantum numbers a — [k', A, fi], 

(r\a) = ^u K , iX (r)Y x>fi (r), (76) 

where we have separated the radial part from the angular part. In the asymp- 
totic limit, we may substitute for the approximate form of the radial part, as 
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derived in Eq. ( 59 1 , 



(f\a) ~ iV K /-sin 



1ii(2kV) Att 

KClQ I 



(77) 



where $ K '.\(r) is the radial-dependent and A-dependent phase-shift, and where 
the coefficient for energy re-normalization is N K i — [2m/(7r 2 ?i 2 e a )] 1 / 4 = 



\J 2m/(irh 

The asymptotic part of the bound states can be neglected, because they decay 
exponentially with the distance from the nucleus and can not affect the outgo- 
ing photoelectron. This allows us to aggregate the sum over the bound state 
spectrum with the contribution of the continuous spectrum in a single integral 
from minus infinity to plus infinity. In this way, the smoothed "bound part" is 
still exponentially damped and it will not contribute, but the integral is easier 
to handle. 



The principal- value part of the wave packet, Eq. (75 1, is evaluated by rewrit 



ing the sine-function in Eq. (77) as outgoing and incoming waves, 



{out/in) 
1 k',X 



(r) 



exp[±i/e'r]. It is then possible to integrate analytically for the principal- value 
using Cauchy's integral theorem, as sketched in Fig. [TTj 




Figure 11: Integration paths in the complex plane for obtaining the principal 
value P = I\ +I2 of outgoing waves: exp[iy/2er] / (e — e K ). The outer integration 
path vanishes, Iq — > 0, as the electron enters the asymptotic region, r — > 00. The 
semicircle over the pole I1/2 becomes — iir times the corresponding residue as the 
arc is reduced. Using Cauchy's integral theorem, we write: I0+I1+I1/2+I2 = 0, 
which implies that the principal value is: P w —I1/2 = ift exp[i-\/2e K r] . A similar 
calculation can be made for the incoming waves, but here the integration paths 
go into the negative imaginary axis which results in an additional minus sign. 
For simplicity, we have here omitted the phase-shifts, & K >,\(r) and the amplitude 
prefactors, since they present no further complication. Atomic units were used 
for compactness. 



The result from this calculation is that the imaginary part of the wave packet 
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takes a cos-like form at the energy-conserving wave number, k, so that 



7T v— * 



+0) y£2(*)> 



(78) 



where the sum over c is needed if the wave packet populates multiple angular 
channels (with short-hand notation: Y^j, so that [A, fi] — [A c ,/i c ]). Inserting 
these asymptotic expressions into Eq. (751, we find that the photoelectron wave 



packet is, indeed, a sum of outgoing Coulomb waves for each angular channel c 

§W(t,f) ps 



neE ■ 



(79) 



Outgoing Coulomb wave 

where the result of the long-range Coulomb interaction and of the short-range 
interactions are contained in ^^\( r ) for the wave number, k, corresponding 
to energy conservation. This intuitive form of the first-order wave packet is 



valid asymptotically, as can be identified in Fig. 12 where we have numerically 
calculated a representative first-order wave packet for hydrogen. 
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Figure 12: First-order radial wave packet (perturbed wavefunction) calculated 
numerically for Hydrogen at a kinetic energy of 40 eV. The real (red) and imag- 
inary (blue) part oscillate as an outgoing wave, as expected from the asymptotic 
approximation. 



Our result for the asymptotic form of the wave packet, Eq. ( 79 ), is quite general 
and it can be applied to many atomic systems by choosing an appropriate phase, 
^k.aM, with different asymptotic phase-shifts, r] Kt \, but also ionic charges, Z. 
Close to the core, the imaginary part of the wave packet will go to zero (it 
will not diverge like the irregular solution [75] )■ The polarization of the bound 
states is clearly visible in Fig. [12] as a "bump" in the imaginary part. Using 
the asymptotic approximation, we neglect this short-range feature, which turns 
out to be an excellent approximation for laser-assisted photoionization, as we 
shall discuss in Sec. [5] The asymptotic form of the wave packet [90] can be also 
expressed in terms of the Coulomb Green's function [HU [H2] . In the special case 
of an hydrogenic system, exact calculations can be performed by using compact 
representations of the Coulomb Greens function [US] [M] . Furthermore, the wave 
packet can be related to the so-called perturbed wavefunction, which satisfies a 
corresponding inhomogeneous, time-independent differential equation [54] 157] . 
In this sense, the wave packet that we have derived using time-dependent per- 
turbation theory is a time-independent quantity, which plays a key role also in 
time-independent perturbative methods! 



4.4 Generalization of first-order photoelectron wave pack- 
ets 

With the help of an integration performed in the complex energy-plane, we have 
just shown that a monochromatic ionizing XUV field, with frequency O, leads to 
an outgoing photoelectron wave packet: j^ 1 ) (t; fl))— (^^(O^expf— iu)(Q)t], 
which oscillates at the characteristic frequency of energy conservation, w(f2) = 
LOi + Q, but with no other time-dependence. Likewise, the effective wave number 
of the wave packet varies with the photon energy so that the phase of the wave 
packet (in angular channel A) is: 



arg 



(80) 



where </>n is the phase of the XUV field, and n(Q) is the wave number upon 
absorption of one such photon with frequency f2. Any attosecond XUV pulse 
can be written as a linear superposition of monochromatic waves, and the cor- 
responding photoelectron wave packets can be written in terms of their time- 
independent (spectral) components, 



1 

2^ 
1 

2^ 



dQ |*W(t;0)) 

dfl exp[-iu(fl)t] . 



(81) 



Time— independent 



Temporal phase 



where ^^(Ojfi) = ^^(il). Note that the superposition presented here is 



different from the one in Eq. ( 45 ) , because the latter includes a superposition of 
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both ingoing and outgoing waves, while Eq. (81 ) only contains outgoing waves 



This is due to the fact that we have already imposed a time-boundary condition, 



namely that the laser field was not on at t — > — oo. As a consequence, Eq. (81 1 
describes the physical process at all times and for any pulse shape, without the 
need for time- dependent complex amplitudes to describe the population of the 



states, as would be required with Eq. (45). The wave packet properties of this 



superposition can be analysed using the SPA, Eq. (31 ), because, in the limit of 
time going to infinity, the oscillations of the phase factor will be fast compared 
to all pre-factors [79]. The SPA is applied as 

* (1) (r,f) = — [ dQ a(fi) exp[if(Sl;r,t)] 
2tt J 

« ^ a(n.) «p[i/(n.)], (82) 

where /(Q;r, t) = /(Q) is the phase function; and a(fi) is a pre-factor, which is 
a slowly varying function of f2. It contains both the spectral envelope function 



of the attosecond pulse and the dipole matrix element. As given in Eq. (80), 
the phase factor for the wave packet (in angular channel A) is 

/(fi; r, t) = n + «(fi)r + 9> K( a),\{r) ~ u(n)t, (83) 

where k and uj are determined by energy conservation from an initial bound state 
after absorbing one photon from the O-field. The stationary phase equation, 



Eq. (32), yields the XUV frequencies, fi s , that give the dominant contribution 



to the electron wave packet at a given space-time position, (r,t) [95]. Using 



Eq. ( 83 ) , the solution to the stationary-phase equation gives the time when an 



electron with wave number k arrives at r 

m = = °> (84) 

TGD T/ree T X +T LR (r) =1 

where we have identified the group delay of the attosecond pulse, tqd(CI); a 
free-particle trajectory r/ ree = v c r; the short-range Wigner delay, t\\ the long- 
range Coulomb delay TXij(r); and the actual time itself, t ^> tqd- This implies 
that the electron with wave vector k will arrive at r at 

t = T GD +T free +T LR + Tx, (85) 

which clearly depends on when the appropriate frequency ionized the atom, tgd] 
and the time it takes for the electron to propagate to the detector, Tf ree + tlr + 
t\. If the electron followed a free-particle trajectory the delay would be Tf ree , 
but due to the Coulomb interaction the electron is delayed by an additional 
amount TLji(r)+T\, corresponding to the long-range and short range corrections 
to the timing of the trajectory. 
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The "delay" from the long-range interaction is difficult to interpret, as we 



already discussed in Sec. 3.3 but we may note that it has a negative sign 
asymptotically 



T LR (r) = - — {1 - ln 2/sr } « -ft^^ 



(86) 



which implies that the electron arrives faster than if it was propagating freely! 



The physical reason is that, as we discussed in relation to Eq. (56), attractive 



potentials lead to a greater local velocity of the electrons. Also the short-range 
interaction may have a negative delay, which classically can be interpreted as 
the electron being somewhat advanced, i.e. starting at a small radial distance, 
a e ff = —t\/v c > 0, outside the origin. In the case of the Coulomb potential, the 
Wigner delays are, however, positive, and such simple arguments do not apply. 
In the asymptotic limit these corrections are very small compared to the free 
travel time and can be neglected. This implies that the dominant contribution 
to the electron wave packet at a given point in space-time arises from 

the part of the wave packet that fulfils the simple relation: v c t sa r. We then 
obtain that the probability density of the photoelectron in an integrated angular 
channel is 

p(i,r) = |* (1) M)| 2 ~k| 2 ^, (87) 

where a K is the complex amplitude of the state with wave number k and the 
time-dependence can be understood as a radial quantum diffusion of the wave 
packet, in good agreement with the simple classical analysis presented above 



in Eq. ( 72 ) . This mapping of complex amplitudes into space and time is the 



basic principle of the TOF detection scheme, which is used in many attosecond 
experiments. In Fig. [13] we sketch the rough dynamics of photoionization. 
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Figure 13: Photoionization in time and space: (a) Initially, at the ionization 
event t = 0, the wave packet is located at the atom. As the electron leaves the 
atom, the slow part of the wave packet is stretched far way from the fast part 
due to quantum diffusion, (b) The momentum distribution, |afc| 2 , of the wave 
packet is naturally accessible using a time of flight (TOF) detector, where the 
probability distribution of fast electrons are detected first, and the slow electrons 
arrive later. To make this example more interesting, the attosecond pulse has 
a "cat" encoded in its momentum distribution, this "cat" eventually appears 
in the space-time distribution. (A more conservative momentum distribution 
would correspond to the odd harmonics of HHG.) 

The photoelectron momentum-to-space mapping, described above, is similar 
to the frequency-to-time mapping for super-chirped light pulses, which was dis- 
cussed in Sec. |2. 2.3) For photoelectrons, the mapping occurs automatically after 
long enough times, and there is no need for a dispersive material. Classically, 
the electron mapping is obvious since electrons in a classical distribution are 
moving at different speeds. If we know that they all started at a distinct point 
in space-time, then we know that the fast electrons will reach the detector (at 
r) first. 

In the case of an APT consisting of high-order harmonics, the photoelectrons 
from the highest harmonic will arrive first, followed by an empty region and then 
the next highest harmonic, and so on. The fact that they were ionized at several 
short instances apart, namely at each peak of the APT, can not be observed 
directly in TOF spectrometer due to the strong quantum diffusion taking place 
in the flight tube. In some sense, the signature of the APT or SAP is seen 
already in the energy and momentum distributions because the electrons do 
arrive in bursts for each harmonic or in a continuous bunch |96j . but in order 
to access the detailed attosecond time scale, one has to probe the photoelectron 
wave-packets with a coherent IR laser field. Next, we shall describe different 
implementations of this concept. 

5 Theory of XUV photoionization 
in the presence of an IR field 

As shown in the preceding sections, attosecond delays in photoionization can be 
deduced from the energy dependence of the phases of the transition amplitudes 
associated with the process. In the experiments reported so far, an auxiliary 
IR laser field is used as a "clock" to measure the delays. Then the analysis is 
somewhat more complicated by the presence of this additional dressing field. 
Indeed, in the presence of an IR laser field, the photoelectron released upon ab- 
sorption of one XUV photon can absorb and emit IR photons, thereby changing 
its energy. Provided that the IR field is weak, the energy of the electron will 
change as the result of the exchange of only one IR photon at most, so that 
the overall process amounts to a two-photon transition. In this section, we give 
a theoretical background relevant to treat this class of laser-assisted photoion- 
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ization, with special emphasis on the contribution from stimulated continuum 
continuum transitions using second-order time-dependent perturbation theory 
[82] . These transitions are associated to the so-called Above-Threshold Ioniza- 
tion (ATI) processes, because the electron is already free when it exchanges the 
second photon. 

The time-dependent dipole interaction with the fields, linearly polarized 
along the z-axis, is written as V(z, t) = Vj(z,t) + Vu(z,t), where we have 



separated the dipole interaction with the XUV field, V[(z, t) in Eq. (66), from 
the subsequent dipole interaction with the laser field, 

Vu(z,t) = ^ (M u ezE u (Sl u )e~ milt jlix. (88) 

With this notation, it is implied that w fli > I p ^> flu ~ u, where the former 
and latter frequencies correspond to the central frequency of the XUV and laser 
field, respectively. We may include a broad bandwidth, Ailj and on both 

fields, provided that the bandwidths do not overlap. The second-order complex 
amplitude for the time-ordered interaction with Vi(t) and Vu(t) is: 



af\t) = 1/ dt' T(f(t')\V n (z,f)\a{t'))a^\t') 



oo 

I 



1/ dt' (f(t')\Vn{z,t')\^(t')), (89) 



where we have identified the first-order wave packet from Eq. (81). It is of 



interest to compare the second-order complex amplitude, Eq. (89), with the first 



order one, as given in Eq. ( 65 1 . In the latter, the electron makes a transition from 
the initial bound state into the field-free continuum, creating a first-order wave 
packet, \^/^(t)). In the second order amplitude, this first-order wave packet 
can be seen as the initial state for the second interaction, and the transition 
towards the final state will occur towards a different continuum state. We can 
write the second-order wave packet formally as 

\¥ 2 Ht)) ee |>£ 2 )(*) \a(t)), (90) 



in analogy with Eq. (64) and our task is now to determine the second-order 



complex amplitudes. After the interactions are over, the S-matrix for the two- 
photon process, Sf' = lim^oo (t), can be rewritten as: 

S% = Jr J dSlj (/H^ 1 )^)) J dn n Eu(Slu) 

/oo 
dt exp[i(wf — flu — u> c (fij))£], (91) 
-oo 

s / 

2it5(uj f —LUi—Q,j—Qjj) 
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where we have expanded the time-dependent quantities as Fourier integrals over 



frequency, using Eq. (81 1 and (66 1, and where we have changed the order of in- 
tegration so that all time-dependence is trapped in a simple inner integral. This 
time-integral is identified as a delta function that enforces energy-conservation 
of the two-photon process, and it is now used to eliminate the spectral integral 
over flu so that the final result becomes: 



fh 



e 1 
ih 2tt 



dSlj </|z|¥«(fi 7 )) EniQ'j 



in 



(92) 



where Vl' n — u)f — Ui — O/. This corresponds to the well-known S-matrix for 
two-photons, here derived assuming a general bandwidth on both fields. 

The two-photon S-matrix, Eq. ( 92 ) , has an intuitive time-independent inter- 



pretation: The integral over fij can be seen as an integral over all the different 
quantum paths that can lead to the final state |/). In order to reach this spe- 
cific final state, the photon energies of the XUV and IR fields must satisfy a 
global energy conservation requirement. If both interactions have a broad band- 
width then there is also a broad range of photon-pairs available. This results in 
a"blurring" effect that arises from the convoluted quantum paths. If one of the 
fields, say the probe, is quasi-monochromatic: fljj = ±uj, which implies that 
Ail j j -c M, the effects of this convolution are negligible. Specializing to the 
dominant contributions associated to the transitions taking place when an XUV 
photon is absorbed first, only two distinct quantum paths are possible: 



'/A 

? (e) 
'//< 



= </H* (1) (0<)) 



ih 
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ih 



(93) 



where E^ = Eu{ui) and E^ = Eji{—uS) = Ejjiui), corresponding to absorp- 
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tion (a) and emission (e) of a laser photon from the fields as depicted in Fig. 
Note that the intermediate XUV photon energy is different in the two paths to 
ensure global energy conservation. 
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Figure 14: Photon- diagrams for laser-assisted photoionization: Single XUV- 
photon absorption contributes with a direct path (d) to the final state. The 
dominant two-photon processes involve the absorption or emission ol a laser 
photon indicated by paths (a) and (e). Higher-order processes involve the ex- 
change ol more laser photons. 

Here, the photoelectron undergoes ATI transitions that require a precise rep- 
resentation of the continuum states, which makes ab-initio computations a 
formidable task for multi-electron atoms or molecules. However, as we shall 
show next, it is feasible to derive a convenient approximation method designed to 
evaluate the phase-shift induced by the continuum-continuum transition stage 
of the process. An accurate determination of these phases are, in fact, essential 
to interpret correctly the delays as measured in recent attosecond time-delay 
experiments. 



5.1 Asymptotic approximation for ATI transition ampli- 
tudes 

In order to obtain an estimate for the phase of the two-photon matrix element, 
we will rely on an approximation, which utilizes the asymptotic forms of the 
final state and the first-order wave packet [JSJ [M] . Both these wave functions 
are of continuum character, with positive energies, and their phase-shifts can be 
derived from their asymptotic behaviour in the limit of a large radial distance, r. 
Atomic units are used in this section for compactness: h = m = e = 1/ (47reo) = 
1. Specializing to the case ol a hydrogenic system with nuclear charge Z, the 
approximate form of the two-photon matrix element for a given photon-pair, 
becomes 

x ^rH)V^ (fe) y iiAf (fc) x Y,( Y L,M\Yi,o\Yx,p)( Y *J Y i,o\Ye, mi ) 

L,M A,/i 
/•oo 

xN k N K drsin[kr + $k,i(r)] r exp[i(nr + $ K , A (r))], (94) 



Asymptotic radial functions 



where we have performed a partial-wave exp ansion with the final state having 



an asymptotic momentum, k, using Eq. (69); and where we have replaced the 
exact form of the radial wave functions by their asymptotic forms. The radial 
functions are written with the short-hand notation for the Coulomb phase-shift, 

. . Zhi(2kr) ttL ,„ , 

**,tW - ^ + Vk,L - -y, (95) 



as defined in Eq. ( 58 >, where r)k,L = o~k,L+$k,L- Here <ik,L — arg{T[X+l— iZ/k]} 
is the pure Coulomb phase with T[z] being the complex gamma function; and 
5k, l is a phase-shift due to an additional short-range interaction. The wave 
numbers of the first-order wave packet, k, and of the final state, fc, satisfy 
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energy conservation after one XUV photon: e K = +Q, and after the exchange 
of one IR photon: efe = e K ± w, respectively. Following Ref. [SI], we evaluate 
the radial integral and obtain: 

% - ~>) 3/ ~ 



E n N k N K 




-iZ(l/n- 
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-i/fe) 


TtZ (\ 




2~ \k 





X 



^2 eX ? 

J] ]T ( Y L, mi \Y lfi \Y Kmi ) 

x {Yx^Y lfi \Y h . m ^{R K MRn^M~ X ^ x{K) - (96) 

Reading this equation from below, the XUV dipolc transition from the initial 
state, \i) = \rii,£i,mi), to the intermediate state, \k), is found on line five. 
The fourth line contains the trivial angular part of the second dipole transition 
from |k) to \k), while the radial part of this continuum-continuum transition 
is divided into two factors on line two and three. On the third line, the pre- 
exponential factor: l/\k — k\ 2 , is characteristic for free-free transitions (59]- If 
the IR frequency is relatively small as compared to the kinetic energy of the 
photoelectron, namely in the soft-photon limit, we have k f=a k — ui/k, so that 
l/\k — k\ 2 ss k 2 /uj 2 , which makes clear that the transition amplitude increases 
for high electron momenta, as well as in the soft-photon limit, lu — > 0. 

Because we consider transitions between Coulomb states, the transition am- 
plitude contains also a real exponential that depends on the nuclear charge, Z, 
together with k and k. More interestingly, the second line is complex and it 
depends on the same three quantities, thus, introducing a new phase-shift into 
the matrix element. Note that these correction factors are universal since they 
depend neither on the angular momentum nor on the short-range atomic in- 
teractions. Finally, the first line contains some pre-factors, such as the electric 
field amplitudes and normalization constants of the wavefunctions associated to 
the relevant electronic states. 

To address the question of the phase of the two-photon matrix element, 
one notices that besides a trivial contribution from the spherical harmonic in 
the final state, Yj,, mi (fc), it contains only phase-shifts that are governed by 
the angular momentum A of the intermediate state, i.e. a state that can be 
reached via smg/e-photon ionization. More precisely, for a given transition chan- 
nel (characterized by the angular momenta of the intermediate and final state: 
£i — > A — > L), the phase of the matrix element reduces to: 

arg[5 (L,A, mi)] = 7r + arg[y i , mi (fc)]+<fo + ^ 

- -y +V\(K) + cf>cc(k,K), (97) 
where </>n and <f>^ are the phases of the XUV field, f2, and of the IR laser, ui, 
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respectively, and where the continuum- continuum phase can be approximated 
by 

r {2k) iZ ^ r[2 + iz(i/ K - i/k)] \ 

E arg \ (2k) iZ / k (k- fc)^(i/ K -i/fe) J ' 



(98) 



corresponding to the argument of line 2 in Eq. (96). The superscript, (P), 
indicates that this result is obtained when using long-range phases of the type 
in Eq. ( 95 1 . The continuum-continuum phases for absorption and stimulated 



emission of one IR laser photon are shown in Fig. |15| 
0.4 
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Figure 15: Continuum-continuum phases: ATI phase-shifts for absorption (a) 
of one 800 nm laser photon; and for stimulated emission (e) of one laser photon, 
in a Coulomb potential with Z = 1. The asymptotic approximation (P) provides 
the correct qualitative behaviour, while the long-range amplitude corrections 
(P+A and P+A') lead to quantitative agreement, at high enough energy, with 
the exact calculations in hydrogen (+ and x). 

We stress that the final state scattering phase, ??i(fc), cancels out and that it 
enters neither in Eq. (96) nor in Eq. (97). 



5.1.1 Long-range amplitude corrections 

So far, we have established that the "long-range phases", i.e. the logarithmic 



divergence of the Coulomb phases in Eq. ( 95 ) , is important for the phase of ATI 
transitions in an ionic potential of charge Z. In order to improve our asymp- 
totic approximation, we must include also long-range amplitude-variations. We 
do this by returning the normalization constants Nk and N K back inside the 
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radial integral in Eq. (94), and write them as functions Nk(r) and N K (r) for 
the radial amplitude variations |54] . In analogy with the local momentum from 
WKB theory, as stated in Eq. ( 50 ) , these effects can then be approximated by 



expanding the normalization factors, so that 



N k (r)N K (r) = 



ir 2 kn 



1- - 



1 



Zi 



(99) 



Correction term 

to first-order in the Coulomb potential. It is also possible to recover a similar 
asymptotic form of the Coulomb functions by using the iterative formulas in 
Ref. |80j . i.e. without the need for a semi-classical approximation. Including 
this correction term, the continuum-continuum phase becomes: 

((2kY z / k T[2 + iZ(l/K-l/k)] + 7(fc,/c)' 
arg \ {2k) iZ / k " (k - jfc)iZ(i/«-i/fc) 

where the additional term is 

(n-k)(n 2 + k 2 ) 



(100) 



j(k, k) = iZ 



2n 2 k 2 



T[l + iZ(l/K-l/k)], 



(101) 



with Eq. ( 100 ) corresponding to the long-range phase and amplitude data (P+A) 
in Fig. |15| We note that this amplitude correction leads to an excellent agree- 
ment with the exact calculation at high energies. The Taylor expansion in 



Eq. (99) is, however, problematic because the Coulomb potential will dominate 
over the kinetic energy close to the core, Z/r > et as r -> 0, which implies 
that the expansion is invalid for r < Z/eu- This leads to an artificial divergence 
of the electron wavefunction close to the core, which is physically unreasonable 
because the wavefunction should vanish at r — 0. 

Because an accurate analytical expression for the continuum-continuum 
phase can be a valuable tool in the analysis of experimental data, we may try 
to avoid the region close to the core, by an ad-hoc mathematical trick, namely 
by changing the starting point of the radial integral. Indeed, the result pre- 
sented as (P+A') in Fig. (15) is slightly better, where the gamma function in 
Eq. ( 101 1 was replaced by an incomplete gamma function represented by an in- 
tegral with an imaginary starting point, ro = iZ(l/ ' k 2 + l/fc 2 )/4. (We note that 
this "regularization procedure" is slightly more systematic than the procedure 
presented in Ref. [54j . but the best matching point was found by trial-and- 
error.) Strictly speaking, at very low kinetic energy, typically below 10 eV, the 
asymptotic approximation starts to break down. Physically, this is due to an in- 
creasing influence of the Coulomb potential for slow electrons. Mathematically, 
this is expected since the asymptotic expressions for Coulomb function have kr 
as arguments, which implies that a small k will lead to a "good" asymptotic 
wavefunction only at a large distance r from the core |80j . 
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5.2 Extracting time-delay information from laser-assisted 
photoionization signals 

Having established the asymptotic approximation for the complex amplitudes of 
ATI processes, we now turn to the probability for the emission of a photoelectron 
with energy = + il, as depicted in Fig. [14] The probability is given by the 
square of the sum of the amplitudes: 

P- k « \S d + S a + S e \ 2 

= |5 d | 2 + |5 Q | 2 + |5 e | 2 

+ 2n{S* d {S a + S e ) + S* a S e }, (102) 
where d, a and e label the paths associated to the lowest-order processes: (d), 



(a) and (e) in Fig. 14 The total probability depends on the relative phase 
of all individual quantum paths and the maximal probability for photoemission 
occurs when all paths are in phase, argfSy = arg[S' a ] = arg[S' e ]. In experiments, 
the phase of the two-photon amplitudes labelled (a) and (e) can be controlled 
by changing the sub-cycle delay, r, between the probe field and the attosecond 
pulses. More precisely, one controls the relative phase of the IR field, <j>± u = 



±ojt in Eq. (97), with respect to the group delay of the XUV pulse. The 



probe-phase dependence is S a cx E u oc exp[zwr] and S e cx £7* cx exp[— iur] 



This implies that the cross-terms in Eq. (102) vary differently as a function 
of t: (d)-(a) and (d)-(e) are modulated with periodicity ujt associated to the 
exchange of only one IR photon [97] : while the cross-term (a)-(e) is modulated 
with periodicity 2ojt due to the two IR photons involved [32j [33] . 

The latter transitions are directly involved in the RABITT scheme, since 
path (d) does not contribute. This is because only odd XUV harmonics are 
used to photoionize the atom, as was indicated in Fig. [5] so that the signal 
from cross-term (a)-(e) occurs at zero-background. Let us emphasize that in 
the case of a SAP of XUV radiation, the delay-dependent modulation of the 
photoionization signal at the central energy O , will also be given by the cross- 
term (a)-(e) alone, because the other two cross-terms cancel with each other 



[53] . All the cross-terms in Eq. (102) together form the on-set of streaking, i.e. 
a small constructive or destructive interference at the high or low momentum 
respectively, so that the average momentum of the photoelectron is deflected by 
different amounts depending on r. Provided that the soft-photon approximation 
is valid, the displacement of this entire structure for a laser-assisted SAP is 
exactly equal to that of APT [S3]. Therefore, the maximal probability for 
photoemission occurs when the amplitudes associated to paths (a) and (e) are 
in phase, arg[S' a ] = arg[S' e ], which is equivalent to 

<p n< +UJT + 7?a(k<) + 4>cc(k, K<) = 

(f> n> -OJT + TJxiKy) + </> cc (fc,K>) (103) 

where we have used the explicit phases of the relevant two-photon matrix ele- 



ments in Eq. ( 97 1 , assuming one dominant intermediate angular channel, with 
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angular momentum A; and momenta k < and n > corresponding to absorption of 
photon il < and £!>, respectively. The solution to Eq. (103) is: 



2lu 

t?a(k>) - 7?a(k<) 4>cc(k,K>) - (j>cc(k, K<) 
2w ~ 2w 



(104) 



where we observe that the probe-delay, r, that maximizes the yield is a sum of 
three delays: 



tqd: the group delay of the XUV field is when the attosecond pulse arrived 
on target, i.e. when the fields J7< and J7> added constructively at the 
atom. 

t\: the Wigner delay is the "delay" in single-photon ionization, i.e. an 
asymptotic temporal-shift of the photoelectron wave packet. Strictly speak- 
ing, this delay is interesting, i.e. giving information on the electron dy- 
namics, only in difference with a Coulomb reference, e.g. hydrogen with 
r] x (k) = <J\(k) = arg{T[l + A — i/n]}. In this way, we may say that the 

( H ) 

photoelectron is delayed by t\ — , as compared to hydrogen. 

t cc : the continuum-continuum delay, i.e. a measurement-induced delay 
due the electron being probed by an IR laser field in a long-range potential 
with a Coulomb tail of charge Z. This delay can be traced back to the 
phase-shifts of the ATI matrix elements. 



We stress that the delays presented in Eq. ( 104 1 are calculated from the finite- 
difference approximations to the actual derivatives, tqd = d4>n/dQ and t\ — 
9T) K> x/dCl. For these approximations to be valid, we must require that the spec- 
tral phases vary slowly, e.g. that |A0n| = \4>n> — 4>n < \ *C 2-7T, corresponding to 
a small phase variation over two laser photons, Afl = 2u. Using the framework 
presented in Sec. |2.1.1| one can easily show that these slowly varying phases 
imply that the frequency components of the attosecond pulse must be confined 
to within a fraction of the laser period. In this sense, "slow reactions" , e.g. reso- 
nances, which induce a dramatic phase jump in the spectral domain are difficult 
to study using the conventional RABITT scheme. In the case of streaking, the 
delays arise also as finite-difference approximations to the actual delays, but 
here the smallest energy difference is instead one laser photon, AO = u> [54] . 
The atomic delays, Tg = t\ + t cc , are shown in Fig. |16| corresponding to the 
case of hydrogen, with t\ for a few different values of the angular momentum, 
A. 
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Figure 16: Atomic delays in laser-assisted ionization: The Wigner delays for 
hydrogen with a photoelectron of s, p, d or /-character [red, upper curves] plus 
the universal continuum-continuum delay for the laser-probing process with 
Z = 1 and Tiuo — 1.55 eV [blue, dashed curve], yields the total delay in laser- 
assisted photoionization [black, lower curves]. 

In this way, it is possible to "measure" the group delay of the attosecond 
pulse provided that the atomic delay: tq = t\+t cc , can be calculated accurately 
and that it is subtracted from the experimental data. The first observation 
of attosecond pulse structures in 2001 was made in this way by subtracting 
the influence of the "atomic phase" as calculated by means of complex two- 
photon matrix elements [31] . We stress that the exact temporal structure of the 
attosecond pulse is not measured directly, neither with the RABITT nor with 
the attosecond streak-camera technique, but that the influence of the atomic 
potential tends to decrease at high kinetic energy. The question of the precise 
determination of these "atomic delays" (also called "streaking delays") and of 
their importance in RABITT and streaking measurements has motivated a large 
number of theoretical studies, see Ref. [HI GU QUI [Ml HQ2J [Ml [Ml [Ml 1M] ■ 
Much of this recent work has focused on the attosecond streak-camera technique, 
where typically non-perturbative approaches are employed to account for high- 
order photon processes with the probe field. In the next subsection, we will also 
comment on the influence of such an intense probe field with a view-point from 
perturbation theory using the asymptotic approximation. 

Alternatively, the Wigner-like delay, t\, can be determined if tqd is known 
from independent measurements and if r cc can be calculated [461 U)T| . However, 
since the precise characterization of the attosecond pulse requires the previ- 
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ous knowledge of t\, it would seem that we have run into an uncomfortable 
circular argument. Instead, if the "same" attosecond pulse is used to photoion- 
ize simultaneously different atomic states, information about the differences in 
atomic delays, Tg = t\ + t cc , can be obtained directly in the experiment, with- 
out any need to know the temporal structure of the attosecond pulse [HJ US] . 



Then, by simply subtracting the universal r cc as given by Eq. (104), the T\- 



difference between the two photoclcctrons is obtained. A similar separation of 



the time-delays as shown in Eq. ( 104 ) was reported by the group of Burgdorfer 
for streaking of SAP by solving numerically the TDSE [104. 105 . Interestingly, 
they also showed that r cc [there called: "the Coulomb-laser coupling"] could 
be calculated from a classical ensemble of electrons interacting simultaneously 
with both the laser field and the Coulomb potential. Prior to these successful 
demonstrations of the seperability of the two contributions, the influence of the 
short-range scattering phase-shifts was shown by Yakovlev and coworkers [99] 
and the importance of the long-range Coulomb tail was discussed by Zhang and 
Thumm [TU2] , 

5.3 Towards higher probe intensities 

So far the discussion has been conducted within a perturbative framework, by 
retaining the lowest-order contributions from the IR probe field. A natural task 
is then to investigate the consequences of an increased probe-field intensity. 
Numerical evaluations of the TDSE show that a strong probe field will alter 
the temporal information [58] . A modest increase of the probe intensity will, 
however, not alter the observed delay. This stability of the delay at low- intensity 
probe fields was first pointed out by Zhang and Thumm [100 within the Eikonal- 
Volkov Approximation [53] by considering a weak probe field. 

Using perturbation theory for the interaction with the probe field makes it 
difficult to draw conclusions about the strong-probe field interaction, but it is 



possible to apply the asymptotic approximation, Eq. (94 1, also for processes in- 
volving additional probe photons. In Fig. |14[ the complex amplitudes involving 
two probe photons are labelled: (d ae ), (d ea ), (aa) and (ee). Here, (d ae ) and (d eo ) 
can be regarded as corrections to the direct path (d), because they exchange two 
probe photons so that the total energy is unchanged. Any such "probe-photon 
loops" , indicated by grey rings in Fig. |14[ involve two continuum-continuum 
transitions: k — > k — > k. Therefore, they acquire two continuum-continuum 
phases: </> cc (k, k) and </> cc (fc, k) = — </> cc (k, fc), which compensate each other. For 
instance, the phase of (d ae ) and (d ea ) is equal to that of the direct path (d) [ex- 
cept for trivial perturbation expansion coefficients]. The two other paths (aa) 
and (ee) yield a new probability modulations at 4wt, due to the four photons 
involved in the (aa)-(ee) cross-term, and they will not disturb the lower-order 
modulations. One can show that the modulation shift of the (aa)-(ee) cross- 
term will be approximately equal to that of the (a)-(e) cross-term, so that no 
fundamentally new information is obtained at this higher-modulation rate, as 
reported in the experimental work by Swoboda and co-workers |107j . In this 
way, including additional probe-photon loops, one can argue for the validity 
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of Eq. (104), also when higher-order contributions become important, at least 
while the probe-interactions remain in the perturbative regime and the bound 
states are not too much dressed. 

In the following section, we will review the state-of-the-art experimental 
efforts to measure these delays in laser-assisted photoionization. As we will 
see, there are so-far few experimental results to compare with the theoretical 
predictions, and in addition, these experimental data points are located in highly 
complex regions, e.g. close to pseudo-resonances in Neon |106j and correlation- 
induced Cooper minima in Argon |108j . where the SAE approximation is likely 
to break down. 



6 Experimental observations of attosecond de- 
lays in photoemission 

The approximations for laser-assisted photoionization described in Sec. [5] hold 
very well for atomic systems, where the photoelectron can be described within 
the SAE approximation. In fact, the accuracy of the asymptotic approxima- 
tion was first bench-marked by comparison with exact calculations in hydrogen 
carried out by R. Tai'eb |46( 154] . The experimental work on photoionization 
time-delays, on the other hand, is carried out on many-electron systems, most 
often noble gas atoms such as Neon and Argon. The question then arises to 
which degree the SAE approximation is valid for these systems, and especially 
so when the photoelectron is released, not from the outer-most orbital, but from 
"inside" the core. Screening effects due to electron correlation can be accounted 
for using many-body perturbation theory (MBPT), such as the random-phase 
approximation (RPA), or more elaborate methods [1011 11051 11061 I108j . but 
we have not considered such corrections here. We predict, however, that in 
the coming years, the interpretation of attosecond experiments will offer a new 
testing ground for many-body calculations. Future experiments will provide 
interesting opportunities to test various theoretical methods, such as MBPT, 
for light-induced electron-electron interactions in connection to measurements 
of phases and delays in complex atomic and molecular systems. At this point, 
we report that a satisfactory agreement between theory and experiment has not 
yet been reached, and that more experiments are needed in "simple" energy 
regions, where perhaps more straight-forward analysis would apply. In the fol- 
lowing, we will proceed with a brief overview of the current experimental efforts. 
We will place special emphasis on the time-delay measurements and phase mea- 
surements using APT, where we have taken active part, but we will also discuss 
the experiments using SAP. We stress that both approaches aim at measuring 
the so-called atomic delays [31], and that the one-photon Wigner delay can be 
accessed only after subtracting the continuum-continuum delay as computed by 
cither the quantum mechanical or classical approaches [J5J 1104) . 
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6.1 Atomic-delay measurements using APT 



A delay of 110 as between the 3p and 3s states in Argon has been measured 
experimentally by Kliinder and co-workers using an APT with a photon energy 
of ~ 35 eV 46] . Theoretically, this delay is identified as a difference of atomic 
phases between two orbitals. These concepts in laser-assisted photoionization 
dates back to the first measurements of attosecond pulses in 2001 using RABITT 
[3TJ [33], but the "atomic delays", by themselves, had never been accessible 
experimentally before due to the simultaneously unknown GD of the attosecond 
pulses. A breakthrough in measuring these atomic delays was possible due to 
the new way that the experiment was performed, where two RABITT scans 
were recorded simultaneously from two different initial orbitals. Interestingly, 
in this setup, the unknown temporal structure of the attosecond pulses can 
be subtracted without ever knowing their exact shape, as can be understood 
by considering Eq. (104 1. While the experiments is distinctly different, the 



method was inspired by the previous delay-measurements using SAP 



A schematic illustration of the experiment is shown in Fig. 17 
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Figure 17: Delay experiment for the 3s and 3p states in Ar using APT: The 
APT is composed of three odd harmonics, 2q + 1, that ionize (pump) electrons 
from either initial state, 3s or 3p, into the continuum. After absorbing a har- 
monic photon, the photoelectron can either absorb or emit a laser probe photon 
in order to reach an even number sideband (SB) state, labelled by 2q. The 
SB probability oscillates with the delay APT-IR delay, r, due to interference 
between the two quantum paths. Information about the attosecond timing is 
found in the relative modulation offset between the same sideband numbers 
from different initial states 

In this work, the high-order harmonic comb was first passed through a thin 
Chrome foil acting as a band-pass to separate four odd harmonics (21-27). In 
this way, only three sidebands, SB: 22, 24, 26, were produced when ionizing ar- 
gon atoms from either orbital. This limiting of the harmonic comb was made to 
prevent different sidebands from different orbitals to overlap in photoelectron ki- 
netic energy. The analysis of the experiment [46] was first carried out within the 
SAE approximation using Hartree-Fock phase-shifts from the litterature [85] , an 
adequate theoretical method for photoelectrons from the outer orbital, 3p. On 
the other hand, photoelectrons from the inner orbital, 3s, couple strongly with 
the outer orbital, and the corresponding single-photon phase-shifts are greatly 
altered. This important effect was identified by Kheifets by including RPA 
effects in the absorption of a single XUV-photon. Interestingly, these photo- 
electrons exhibit a correlation-induced minimum in the one-photon ionization 
step, i.e. a complex kind of Cooper minimum [84,, with a corresponding peak 
of several hundreds of attoseconds in the atomic delay |108j . This is, indeed, a 
challenging spectral region, where several quantum paths can interfere, and the 
results become difficult to interpret theoretically. 

The signal can also be strongly altered by atomic resonances. In the spec- 
tral domain, an APT corresponds to harmonics that can be "aimed" towards 
specific energy regions of interest. The frequency of the harmonics from the 
HHG process can be tuned naturally by changing the fundamental (driving) 
laser pulse frequency, w + Slo, so that the high-order harmonic frequencies in- 
crease or decrease by Sfl = (2q + l)Sw. Swoboda et al. found that the phase of 
the modulation of the lowest-lying sideband in Helium depended critically on 



the frequency of the harmonics |109j . In Fig. 18 an illustration of the principle 
of the experiment is shown. 
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Figure 18: Phase induced by the ls3p state in He using APT: The phase 
in photoionization, 5t, depends on the detuning, Sfl, between the harmonic, 
(2q — 1)lo = 15cj, and the bound state, ls3p. The effect arises from the 7r-jump 
that occurs in the two-photon amplitude, when passing the resonance. The 
higher-lying sideband, SB: 2q + 2, acts as an important "reference clock" in 
this experiment, which is mostly insensitive to small frequency changes of the 
harmonic fields [109] . 

The observed modulation-shifts were attributed to a phase effect occurring over 
the intermediate bound state, ls3p 1 P±, which was "on resonance" with a below- 
threshold harmonic, 157kJ < Ip He \ The modulation-shift was then compared 
to the corresponding case in Argon, where the harmonic instead directly creates 

(Ar) 

a photoelectron, 15Jiuj > Ip , and where the bound-states do not contribute 
considerably. By tuning the harmonic over the sharp resonance, an abrupt ir- 
shift is expected. The experimental 7r-shift is, however, "smoothed" due to the 
finite duration of the laser probe field and of the APT (~ 30 fs), which is much 
shorter than that of the long-lived bound state. In order to fully resolve the 
abrupt rise of the delay due to pure atomic effects, one needs harmonics that 
are more narrow in frequency than the atomic resonance, i.e. an APT which 
is longer that the lifetime of the state. Similar phase-shifts have also been 
found and analysed at short-lived complex continuum resonances in Nitrogen 
molecules |110l llllj . In this case, the life-times are more comparable to the 
duration of the APT, implying that observed modulation-shifts contain some 
information of about the lifetime of the resonances. 
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6.2 Atomic-delay measurements using SAP 



It was for a long time believed that the so-called "attosecond streak-camera" 
|35j , provided a direct link between the time-domain and the energy-domain, as 
often illustrated by measuring the oscillations of few-cycle optical light pulses 
[55] . This stand-point needed to be revised after an experiment by Schultze and 
co-workers, where it was stated that "the zero of time in atomic chronoscopy 
is currently tainted with an error of up to several tens of attoseconds" [45] . 
due to a measured time-delay of ~ 20 as between the streaked photoelectrons 
from the 2p and 2s shells in Neon at an XUV photon energy of 100 eV. The 
key to obtaining this intriguing experimental result was to collect simultane- 
ously two streaked photoelectron spectrograms from two different initial atomic 
states, thus accessing the difference in streaking delay between the two pro- 
cesses. Interestingly, by measuring this difference in timing, the unknown GD 
of the attosecond pump pulse cancels out, because it is the same in both pro- 
cesses, as can be identified using Eq. (104 1. The experimental results motivated 
at great deal of further theoretical work, where it was concluded that the long- 
range Coulomb tail is responsible for a considerable part of the observed delay 
in streaking [100[ I102[ 1104] . In addition to this long-range effect, the laser- 
induced polarization of the atom can lead to a shift of the streaking delay if the 
probe field is too strong (S3J Extensive theoretical work including many 
multi-electron effects have so far only accounted for a fraction of the experi- 
mental time-zero shift [1011 1106] . The theory of the attosecond streak-camera 
has been refined to account for the long-range Coulomb tail by matching the 
asymptotic wavefunction to the appropriate scattering phase [103] . but these 
modifications are controversial as they go against Eq. (104 1, as is also stressed 
in Ref. [105] , More experimental results from other systems than Neon, are 
required to compare with these new theoretical models. 

In conclusion, the attosecond pump-probe schemes using a SAP or APT, 
Fig. [4] (b) and (c) respectively, can be used to obtain information about the 
probe field-dressed system, and not directly the probe-free system. To explore 
properties of the unprobed system using experimental data where the probe is 
present, a theoretically reliable way of accounting for the probe must be found 
[4*51 1104] , By decreasing the intensity of the probe field, the streaking modula- 
tion is reduced and one arrives in a regime where standard perturbation theory 
is applicable [52] • This presents a large advantage for quantitative theoretical 
work, since it can be based on well-established MBPT. 



6.3 Other attosecond time-delay experiments 

Finally, in this subsection we briefly discuss other attosecond delay experiments 
that are not directly related to single-photon ionization by an XUV field. This 
important class of attosecond experiments are conducted in the HHG cell di- 
rectly, and they are sometimes referred to as self-probing systems, see Ref. [112] 
for a recent tutorial on self-probing of molecules with HHG. In these experi- 
ments, the electron wave packet in the continuum is thought of as a probe for 
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the entire system, holding promise of simultaneous temporal and spatial reso- 
lution of molecular dynamics [50] . In this subsection, we will restrict ourselves 
to one of the simplest self-probing systems, namely, atoms undergoing HHG in 
two-colour laser fields consisting of a strong fundamental, u, and a much weaker, 
parallel second harmonic, 2oj. An interesting aspect of this two-colour setup, 
is that the presence of the second harmonic breaks the inversion symmetry of 
the HHG process, which leads to the production of odd and even harmonics 
|113] 1114] 1115] . A phase-difference (or time-delay) in the intensity modulation 
of different even harmonics was observed as a function of the delay between the 
to and 2uj fields. Early experimental and theoretical work suggested that the 
phase of this modulation could be used to characterize the emission times of 
attosecond pulses [3TJ [551 1116] . Experimentally, it was eventually found that 
the phase behaved in an unexpected way close to the harmonic cut-off [117] . 
which, after some more detailed calculations [72], lead to a new interpretation 
of the observed delays with, instead, a connection to the Keldysh tunnelling 
parameter [118] , 

The delay in tunnelling ionization has been studied by the group of Keller 
using angular streaking, often referred to as the "attoclock" |119l 11201 1121] . 
In this setup, a close-to-circularly polarized femtosecond laser field is used to 
tunnel- ionize electrons. Subcycle resolution is obtained by studying angular 
displacements of the final momentum of the photoelectrons. Using this angular 
streaking method a vanishing tunnelling delay has been confirmed. We need 
to stress the differences between the two kinds of streaking: The conventional 
form of streaking requires a SAP and a femtosecond probing laser field with lin- 
ear polarization, while the angular streaking requires only a single femtosecond 
pulse with a tailored polarization. In regular streaking, temporal information 
about laser-assisted, single-XUV-photon ionization is gained; while in the case 
of angular streaking, the temporal aspects of strong-field (multi-photon) tun- 
nelling is studied. The theoretical work presented in Ref. |121) highlights the 
difficulties and prospects of self-probing systems involving electron continuum 
probes initiated through strong-field tunnelling. 

7 Conclusions and outlook 

Attosecond science is attracting a great deal of attention because it promises ac- 
curate control and probing of electron processes in atomic and molecular systems 
in real time. In this tutorial, we have presented a detailed theoretical analysis 
of photoionization of neutral atoms by attosecond XUV pulses and coherent 
XUV high-order harmonics. We have found that, indeed, attosecond pulses can 
be used to initiate photoionization at well-defined times, but to probe such an 
ultra-fast event requires careful considerations about the quantum mechanical 
interactions between the photoelectron and the remaining ion. Ultimately, it 
is commonly believed, that the probing process will be carried out using a sec- 
ond attosecond pulse, but such experiments are demanding due to the need for 
two-photon processes in the XUV range. We refer the reader to Ref. [51] for 



a state-of-the-art experimental scheme using XUV for both pumping and for 
probing the processes on a time-scale bordering the attosecond domain. 

At the present time, the temporal aspects of photoionization are more often 
probed by a phase-locked, IR-laser field. This has the advantage of a much 
stronger interaction with the photoelectron, in fact, here the interaction can be 
so strong that many such laser photons are absorbed. This situation is called 
"streaking" because it can shift the photoelectron momentum distribution to 
higher or lower values depending on the subcycle delay between the attosecond 
pulses and the probe field. This highly non-linear process is difficult to handle 
with high-fidelity in theory, and this is the main reason why we have considered 
instead the interaction with a moderate probe field, which may induce exchange 
of a single IR-laser probe photon at most. 

Such moderate probe fields are the standard choice in the temporal prob- 
ing of APT using the RABITT scheme. We have demonstrated [54] that this 
interferometric interpretation is valid also for SAP in the on-set of streaking. 
Using this "photon-picture", we have identified the phases of the lowest-order 
quantum paths and we have shown that the same temporal information about 
the photoionization process is gained with either pulse structure, SAP or APT. 

Placing an interpretation on the "delay" in laser-assisted photoionization, is 
not as straight-forward as indicated by the early streak-camera formalism |36j . 
This is because the strong-field approximation is not a reliable approximation 
for laser-assisted photoionization. Instead, it turns out that the long-range 
interaction between the photoelectron and the remaining ion leads to a shift of 
the time-delays (100] . which often dominates over the one-photon, Wigner-like 
delays of the photoelectron [99] . Using time-dependent perturbation theory, 
we have shown that these "measurement artifacts" can be separated from the 
Wigner-like delays, as written explicitly in Eq. (104 1, in the case of both APT 
jUj and SAP [104] . Furthermore, we have shown that these temporal artifacts 
can be traced back to phase-shifts in ATI transition-matrix elements between 
continuum states, hence the name "continuum-continuum phases" [54]. We 
have here presented both correction curves and analytical expressions for the 
purpose of correcting for such artifacts in the quest for the true delay in single- 
photon ionization. 

In this tutorial, we have limited ourselves to the SAE approximation, with 
a single photoelectron in a static atomic potential. While this approach is 
certainly valid for hydrogenic systems [SI], it is likely to break down in many- 
electron atoms, when electron-hole interactions are strong and when the electron 
relaxation plays an important role for the dynamics. Here, accurate calculations 
that include multi-electron screening effects of the ion are required. Although 
a few full-scale simulations have been made including many-electron dynamics 
|10R 1106) , a conclusive agreement between experiment and theory is yet to be 
presented. This calls for both more experimental data and detailed theoretical 
work, and then a careful comparison of the two. In this way, the time-delay 
experiments provide an important testing ground for more exact analysis and 
interpretations in attosecond physics and strong-field physics. Only when under- 
standing these detailed ultra-fast phenomena, will "attophysics" truely provide 
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a route to accurate control and probing of electron processes in atomic and 
molecular systems in real time. 
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